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Abstract 


IBM  Reference  No. — TR  22. 1486 

This  report  outlines  results  of  two-dimensional  mathematical  inves¬ 
tigations  initiated  during  the  present  contract  period:  the  insulated 
gate  field-effect  transistor,  the  Schottky  barrier  field-effect  tran¬ 
sistor  and  the  bipolar  transistor.  In  addition,  information  is  presented 
on  the  mathematical  technique  used  to  solve  these  semiconductor  device 
problems.  The  second  investigation  (Schottky  barrier  field-effect 
transistor)  is  complete  and  has  been  published  in  the  technical  literature. 
The  other  two  investigations  (insulated  gate  field-effect  transistor  and 
the  bipolar  transistor)  remain  incomplete,  although  the  material  presented 
here  shows  that  substantial  modification  is  required  in  the  available 
theory  of  these  semiconductor  devices. 

A  discussion  is  presented  on  a  newly  developed  method  for  numer¬ 
ically  solving  the  two-dimensional  ambipolar  diffusion  equations  for 
holes  and  electrons  in  semiconductor  material .  All  numerical  computations 
of  this  type  utilize  underrelaxation  as  a  means  to  attain  computational 
stability;  the  present  discussion  outlines  the  automation  of  this  under¬ 
relaxation  process  whereby  all  device  calculations  can  be  performed  on  a 
"hands-off"  basis.  Another  new  aspect  of  this  computational  method  is  the 
adoption  of  a  formulation  that  provides  rapid  convergence  of  the  cal¬ 
culated  electric  current  in  investigations  of  semiconductor  devices. 
Resulting  from  this  computational  method  is  an  ability  to  calculate  the 
volt-ampere  characteristics  of  semiconductor  devices  with  a  minimum 
number  of  simplifying  assumptions. 

In  two  spatial  dimensions,  calculations  have  been  performed  to 
demonstrate  that  a  soft  current  saturation  is  characteristic  of  JFET 
structures  containing  a  low  conductance  semiconductor  substrate.  The 
calculated  volt-ampere  characteristics  of  this  device  are  presented  in 
conjunction  with  the  mobile  carrier  distribution  that  results  when  the 
structure  is  biased  well  into  electric  current  saturation.  These 
Investigations  include  the  consequences  of  a  velocity  limited  transport 
of  carriers  within  the  source-drain  channel  of  a  JFET.  Resulting  from 
these  calculations  is  a  mathematical  evaluation  of  the  influence  a  low 
conductivity  substrate  has  upon  the  electrical  properties  of  this  semi¬ 
conductor  device;  this  topic  is  of  particular  importance,  since  the 
low  conductivity  substrate  material  is  vulnerable  to  radiation. 


Numerous  solutions  have  been  completed  for  a  boundary  value  problem 
that  closely  approximate  the  two-dimensional  cross  section  of  an  insulated 
gate  field-effect  transistor.  These  solutions  have  established  the 
detailed  steady-state  mechanisms  of  operation  within  this  semiconductor 
device.  For  example,  it  is  shown  that  channel  pinch-off,  as  outlined  in 
elementary  IGFET  theory,  cannot  be  verified  by  a  rigorous  mathematical 
analysis.  An  Important  result  of  this  effort  is  a  detailed  understand¬ 
ing  of  the  physical  mechanisms  producing  electric  current  saturation  in 
an  IGFET,  and  a  recognition  that  the  concept  of  a  well  defined  threshold 
voltage  cannot  be  established  In  the  same  manner  it  has  ^een  previously 
established  from  one-dimensional  IGFET  theory.  These  calculations  also 
provide  insight  into  some  problems  associated  with  our  present  concept  of 
electrical  conduction  in  a  source-drain  channel.  Computational  experiments 
have  been  used  to  demonstrate  basic  differences  between  the  physical  mec¬ 
hanisms  of  IGFET  operation,  as  established  through  a  rigorous  two-dimen¬ 
sional  solution  of  this  problem,  and  as  determined  from  theoretical  concepts 
presently  available  in  the  technical  literature. 

In  the  theory  of  bipolar  transistor  operation,  substantial  importance 
Is  given  to  high  current  hpE  fall -off  and  emitter  current  crowding.  Present 
theoretical  concepts  of  these  observed  characteristics  are  based  upon  one- 
dimensional  and  pseudo  two-dimensional  solutions  of  the  bipolar  transistor; 
there  is  little  (or  no)  verification  of  these  solutions,  using  a  rigorous 
two-dimensional  analysis  of  the  problem.  Froi  work  conducted  under  this 
contract,  it  is  shown  that  substantial  hp£  fall-off  is  attributable  to  two- 
dimensional  mechanisms  involved  in  transistor  operation.  Furthermore,  it 
is  shown  that  one  consequence  of  these  mechanisms  is  that  substantially  less 
emitter  current  crowding  is  observed  than  would  be  predicted  using  available 
bipolar  transistor  theory.  One  consequence  of  this  effort  is  a  realization 
that  we  must  reevaluate  many  mechanisms  that  heretofore  proposed  that  the 
high  current  fall -off  of  current  gain  and  ft  is  a  consequence  of  current 
crowding;  a  detailed  understanding  of  this  hpE  and  ft  problem  is  not 
presently  available. 
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With  partial  support  from  Air  Force  Cambridge  Research  Labor¬ 
atories  under  contract  F  19628-70-C0Q98,  Project  4608,  the  Inter¬ 
national  Business  Machines  Corp.,  Components  Division,  E.  Fishkill, 

N.Y.  has  conductedrmathematical  investigations  directed  toward  a 
study  of  the  influence  of  ionizing  radiation  upon  the  operation  of 
semiconductor  devices.  A  principal  objective  of  this  effort  is  to 
desian  and  develop  computational  techniques  whereby,  in  two  spatial 
dimensions,  the  ambipolar  diffusion  equations  for  holes  and  electrons 
in  a  semiconductor  can  be  solved  for  numerous  different  semiconductor 
devices.  This  effort,  has  been  undertaken  with  the  realization  that  a 
detailed' knowledge  of  semiconductor  device  operation  is  an  important 
step  toward  understanding  the  consequences  of  operating  these  devices 
in  a  radiation  environment. 

^In  conjunction  with  the  development  of  these  computational 
techniques,  problems  of  semiconductor  device  operation  have  been  under¬ 
taken  ( under  this  contract.  These  problems  were  selected  on  the  basis 
of  mutual  interest  to  the  Air  Force  and  to  International  Business 
Machines  Corp.  By  undertaking  these  problems,  new  computational  tech¬ 
niques  have  been  developed  which  greatly  simplify  the  process  of 
mathematically  analyzing  the  operation  of  semiconductor  devices.  In 
addition,  extensive  insight  has  been  gained  into  the  mechanisms  of 
operation  within  semiconductor  devices. 

'  A  substantial  effort  has  been  directed  toward  a  two-dimensional 
mathematical  analysis  of  IGFET  operation.  This  analysis  was  used  as 
a  basic  vehicle  to  refine  the  computational  techniques  developed  under 
this  contract  and,  in  addition,  to  gain  a  detailed  understanding  of 
the  mechanisms  involved  in  the  operation  of  this  semiconductor  device. 

A  principal  goal  for  this  effort  was  to  develop  two-dimensional  com¬ 
putational  methods  whereby,  the  volt-ampere  characteristics  of  this 
structure  can  be  calculated  with  sufficient  accuracy,  and  using  a 
reasonable  amount  of  computer  time. 


In  addition-to^this  investigation  of  IGFET  operafior£,_a  two-dimen¬ 
sional  study  was  made  concerning  the  physical  mechanisms  involved  in 
the  oper»tic^  of  JFET  structures  containing  a  low  conductivity  sub¬ 
strate.  ^This  subject  is  of  importance  to  a  particular  class  of 
integrated  circuits  because  it  provides  substantial,  insight  into  the 
role  played  by  the  conductivity  substrate  material,  which  is  par¬ 
ticularly  vulnerable  to  radiation. 

Another  application  of  these  computational  methods ,-<ieve4eped 
dur4«g-^frirpHrH)d>  is  a  two-dimensional  solution  of  the  bipolar  tran¬ 
sistor.  There  are  two  purpose?  for  undertaking  the  solution  of  this 
problem:  first,  to  test  the  numerous  pseudo  two-dimensional  solutions 
that  are  described  in  the  technical  literature  and,  second,  to  establish 
that  this  extremely  difficult  probl?*m  can  be  solved  In  a  rigorous  fashion 
using  presently  available  high-speed  electronic  computers.  Much  of  the 
radiation  hardening  work  under  way  in  this  country  is  based  upon 
theoretical  calculations  using  pseudo  two-dimensional  techniques,  despite 
the  fact  that  tnese  techniques  have  not  been  subjected  to  a  rigorous 
evaluation. 
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CHAPTER  I 


Two-Dimensional  Analysis  of  JFET  Structures 
Containing  a  Substrate 

1.0  Introduction 

In  a  recent  publication  M.  Reiser1  discussed  the  results  of  a  two- 

dimensional  calculation  of  mechanisms  involved  In  the  operation  of  JFET 

structures  containing  a  low  conductivity  substrate.  His  calculations 

show  that  the  substrate  produces  a  reduction  of  saturation  resistance, 

when  compared  with  the  saturation  resistance  of  devices  containing  an 

2 

Ideally  Insulating  substrate.  This  reduction  of  saturation  resistance 
Is  attributed  to  an  Increased  transverse  component  of  the  source-drain 
channel  electric  field  (with  an  increase  of  source-drain  voltage)  that 
"weakens  the  diffusion  field  below  the  source  and  gate  so  that  more 
carriers  are  pulled  into  the  substrate".  This  concept  is  reinforced 
by  illustrations  of  the  calculated  increase  of  mobile  carriers  within 
the  substrate,  when  this  structure  is  biased  well  Into  electric  current 
saturation. 


The  purpose  of  this  chapter  is  to  present  a  more  detailed  comparison 
between  the  physical  mechanisms  of  operation  within  JFET  structures  con¬ 
taining  an  Ideally  insulating  substrate  and  v- i thin  structures  containing 
a  low  conductivity  semiconductor  substrate.  A  conclusion  derived  from 
this  investigation  Is  that  carrier  accumulation  does,  indeed,  take  place 
within  the  substrate  material  of  this  semiconductor  device,  and  the  mec¬ 
hanist.  producing  accumulation  are  attributable  to  the  field  dependent 
mobility  ef  carriers  within  the  source-drain  channel. 


1.1  Mathematical  Methods 


It  has  been  shown  that  the  hole  and  electron  distributions  in  semi- 

3 

conductor  material  are  described  by  the  equations 


a) 

b) 

c) 

l) 

a) 

f) 


div  grad  ^  (N  «  n  +  p) 

iceQ 

3p  =~  9Dp  9rad  P  '  %P  9rad  * 

3n  =  qDn  grad  n  -  qpnn  grad 

div  iL  =  qR 
P  P 

div  iJ  =  qR„  and 
n  ^  n 

+  3n 


(1) 


1 


by  assuming  no  trapping  mechanisms  within  the  structure  under  consideration. 

*  1  / 
Equation  (la)  is  Poisson's  equation,  and  it  relates  the  divergence  of 

the  electric  field  (E  *  -  grad  ip)  to  the  electrostatic  charge  distribution  1 

arising  from  both  mobile  charge1 carriers  (holes  p  and  electrons  n)  and 

immobile  Ionized  Impurity  atoms  (N)  within  the  semiconductor  lattice. 

1  •  ■  !  * 

Equations  (lb)  and  (1c)  give  the  electric  current  distribution  In  a 
semiconductor  arising  from  the  transport  of  mobile  holes  and  electrons. 

These  equations  express  the  dependency  of  each 'electric  current  component 
(3p  and  $n)  upon  the  concentration  gradients  pf  holes  and  electrons,  the  ! 
mobility  of  these  charge  carriers*  and  the  electrostatic  potential  gradient 
(electric  field)  within  che  semiconductor  material.  1 

i  , 

Equation  (I'd)  anp  (la),  the  continuity  equations  for  holes  and  elec-  .  1 

1 

trons  in  a  semiconductor,  are  based  upon  an  unspecified  mechanism  for 
carrier  generation  and  recombination.  For  simplicity,  the  present 
analysis  is  based  upon  the  assumption  of  an  infinite  minority  carrier 
lifetime  within  the  semiconductor  material.  1  This  simplification  is 
assumed  reasonable  because  the  magnitude  of  recombination  and1  generation  i 
encountered  in  modern  silicon  devices  has  little  influence  upon  the 
Important  mechanisms  contributing  to  field-effect  device  operation. 

Equation  (If).  states  that  the  total  elpctric  current  density  (3y) 
is  a  vector  sum  of  the  electric  current  densities  due  to  holes  (Jp) 
and  electrons  (Jn). 

Finite  difference  methods  are  used  to  solve  this  system  of  equations 
for  the  boundary  value  problems  used  tp  approximate  a  JFET.  In  particular, 
the  current  equations  Qb  and  lc)  are  transformed  into  Mfferential  equa- 

‘  *  I 

tions  for  electrical  flow  lines  within  this  semiconductor  s.-ucture,  and 
these  flow  line  equations  are  solved  (in  conjunction  with  Poisson  ' 
equation)  for  c  large  array  of  spatial  locations.  The  details  of  this 
computational  technique  are  the  topic  of  discussion  in  a  later  chapter. 

Throughout  this  investigation  two  different  JFET  structures  are  con¬ 
sidered:  a  JFET  in  which  the  junction  "pinches-off"  against  an  ideal 
insulating  surface  (Fig.  1),  and  a  JFET  in  which  the  junction  "pinches- 
off"  against  a  substrate  of  low-doped  semiconductor  material  (Fig. 2). 
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The  technique  for  mathematically  modeling  Fig.  1  was  discussed  In 

2  4 

previous  publications  on  this  topic.  *  For  Fig.  2,  a  nodal  array 
was  selected  In  which  the  point  density  is  significantly  increased 
throughout  regions  requiring  detailed  mathematical  Investigation.  In 
particular,  an  Increase  of  point  density  is  used  throughout  the  boun¬ 
dary  between  high-doped  and  low-doped  semiconductor  material,  and 
throughout  that  region  of  the  structure  in  which  the  gate  junction 
"plnches-off"  electric  current  between  the  source  and  drain  contacts. 

In  this  analysis  particular  emphasis  is  placed  upon  the  selection 
of  boundary  conditions  that  do  not  introduce  errors  In  the  calculated 
results.  For  example,  the  boundary  conditions  used  in  this  analysis  ap¬ 
proximate  the  physical  and  electrical  characteristics  of  the  outer 
periphery  of  each  semiconductor  structure,  rather  than  Internal  boundaries 
established  by  Its  physical  or  electrical  properties.  A  homogeneous 
impurity  atom  distribution  is  assumed  within  each  region  of  Figs.  1  and 
2  (gate,  channel,  substrate),  with  an  abrupt  Impurity  transition  at  the 
metallurgical  boundary  of  the  gate  junctions  and  at  the  channel -to-sub- 
strate  boundary  shown  In  Fig.  2.  The  exposed  semiconductor  surfaces  are 
assumed  to  be  an  ideal  Insulator,  i.e.,  no  electric  current  Is  permitted 
normal  to  these  bounding  surfaces.  The  ohmic  contacts  In  each  model  are 
approximated  by  equipotantial  surfaces  that  are  charge  neutral;  these 
•ohmic  contacts  have  been  located  a  sufficient  distance  from  the  active 
regions  of  these  devices  to  have  no  influence  upon  their  electrical 
characteristics. 

1.2  The  High-low  Semiconductor  Junction 

The  channel -to-substrate  boundary  in  F.ig.  2  is  a  high-low  semicon¬ 
ductor  junction;  this  topic  has  been  discussed  extensively  in  the 
technical  literature.^  Because  many  aspects  associated  with  the  operation 
of  this  JFET  structure  (Fig.  2)  involve  electrical  conduction  in  a  high- 
low  junction,  a  brief  review  of  this  subject  is  presented  here.  It  will 
be  shown  that  the  majority  carrier  accumulation  calculated  by  M.  Reiser^ 
is  not  attributable  to  conventional  mechanisms  of  electrical  conduction 
across  the  channel  and  substrate  boundary.  Majority  carrier  accumulation 
in  the  low-doped  region  of  a  high-low  junction  is  a  consequence  of 
conductivity  modulation  due  to  the  accumulation  of  minority  carriers. 


i 

i 

In  the  calculations  of  H.  Reiser,  minority  carriers  were  assumed  to  have 
little  (or  no)  influence  upon  the  mechanisms  of  operation  within  this 
semiconductor  device;  for  this  reason,  the  equation  for  minority  carriers 
was  not  included  in  his  analysis.  Therefore,  conventional  conductivity 
modulation  (and  hence  conventional  majority  carrier  accumulation)  cannot 
take  plac'  in  the  mathematical  model  used  M.  Reiser  in  his  publication 
on  this  subject. 

Figure  3  shows  the  calculated  majority  carrier  distribution  in  an 
n+-n  high-low  junction  at  potential  equilibrium.  From  Fig.  3,  an  abrupt 
transition  between  the  high-doped  and  low-doped  material  is  not  represen¬ 
tative  of  the  mobile  carrier  distribution  at  this  boundary.  An  electro¬ 
static  double-layer  Is  formed  by  a  redistribution  of  majority  carriers, 
relative  to  the  distribution  of  ionized  impurity  atoms.  Carriers  are 
depleted  from  the  high-doped  region  thereby  forming  the  positive  half 
of  the  double-layer;  the  negative  half  is  formed  by  these  same  carriers 
accumulating  on  the  low-doped  side  of  the  high -low  boundary.  Unlike  a 
p-n  junction,  the  electrostatic  charge  produced  by  this  double-layer  is 
not  significantly  influenced  by  the  application  of  an  external  biasing 
voltage.  An  excessively  large  biasing  voltage  (hence,  an  extremely 
large  current  density)  would  be  required  to  significantly  modify  this 
potential  barrier  without  taking  into  consideration  the  mechanism  of 
minority  carrier  conduction. 

Electrical  conduction  across  the  high-low  junction  in  this  JFET 
involves  mechanisms  associated  with  boch  a  forward  biased  and  reverse 
biased  structure.  Specifically,  if  the  external  biasing  voltage  is 
applied  in  a  direction  that  would  tend  to  reduce  the  potential  barrier 
in  a  high-low  junction  (majority  carriers  flow  from  the  high-doped  to  the 
low-doped  material)  the  structure  is  presumed  to  be  forward  biased. 
Similarly,  if  this  applied  biasing  voltage  tends  to  increase  this  poten¬ 
tial  barrier  (majority  carrier  flow  from  the  low-doped  to  the  high-doped 
material)  the  structure  is  presumed  to  be  reverse  biased.  These  two 
different  modes  of  operation  have  an  opposite  influence  upon  the 
majority  carrier  density  in  the  low-doped  region:  forward  bias  produces 
carrier  accumulation  and  reverse  bias  produces  carrier  depletion. 
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Fig.  3.  Impurity  atom  and  mobile  electron  distribution 
in  an  n  -  n  type  of  high-low  junction. 
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In  a  forward  biased  high-low  junction  the  electric  field  distribution 
at  locations  removed  from  the  double-layer  is  determined  by  the  majority 
carrier  distribution  and  the  electric  current;  this  electric  field  distri¬ 
bution  produces  electric  current  continuity,  due  to  majority  carriers, 
throughout  this  entire  semiconductor  structure.  Majority  carrier  electric 
current  continuity  does  not  imply  a  similar  situation  of  electric  current 
continuity  for  minority  carriers,  without  the  intervention  of  diffusion 
mechanisms.  A  detailed  analysis  of  this  structure  shows  that  minority 
carrier  transport  in  the  high-doped  region  is  attributable  almost  entirely 
to  diffusion  from  the  high-low  junction  to  the  high-doped  region  ohmic 
contact.  In  contrast,  minority  carrier  current  in  the  low-doped  region 
is  attributable  to  both  drift  and  diffusion.  Minority  carriers  drift 
toward  the  high-low  junction,  although  the  electric  current  due  to  drift 
far  exceeds  the  electric  current  due  to  diffusion  on  the  h’gh-doped  side; 
hence,  carrier  accumulation  takes  place  in  the  low-doped  material.  This 
process  of  minority  carrier  accumulation  produces  a  diffusion  gradient 
in  a  direction  that  oppose,  the  drift  component  of  minority  carrier 
current,  a*d  thereby  current  continuity  is  obtained  for  minority  carriers. 
This  process  of  minority  carrier  accumulation  will  conductivity  modulate 
the  low-doped  side  and  produce  an  accumulation  of  majority  carriers. 

This  process  of  carrier  accumulation  is  significantly  influenced  by 
the  gradient  of  minority  carriers  oa  the  high-doped  side;  an  increase 
of  gradient  produces  an  increase  of  minority  carrier  current  and,  hence, 
a  decrease  of  minority  carrier  accumulation  in  the  low-doped  material. 

A  consequence  of  this  mechanism  is  a  decrease  of  accumulation  as  the 
gate  junction  "pinches-off"  against'  the  high-low  junction  in  Fig.  2, 
due  to  an  increase  of  minority  carrier  current  produced  by  the  gate  junction 
space-charge  layer. 

Figure  4  illustrates  the  calculated  majority  carrier  accumulation 
within  a  forward  biased  high-low  junction.  Figure  5  shows  the  calculated 
minority  carrier  distribution  producing  this  majority  carrier  accumulation 
(Fig.  4);  in  addition,  Fig.  5  shows  minority  carrier  accumulation  in  the 
high-doped  side  of  this  structure.  Accumulation  of  minority  carriers  in 
the  high-doped  region  is  insufficient  to  produce  a  significant  amount  of 
conductivity  modulation. 


For  a  reverse  biased  high-low  junction,  the  mechanisms  producing 
minority  carrier  depletion  differ  little  from  the  mechanisms  producing 
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Calculated  minority  carrier  distribution  in 
forward  biased  high-lav  junction  (silicon). 
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minority  carrier  accumulation  in  a  forward  biased  stricture.  As  before, 
minority  carrier  current  continuity  must  exist  throughout  the  entire 
high-low  junction.  'Because  the  electric  field  distribution  is  estab-, 
lished  by  majority  carriers,  the  procest  of  drift  alone  cannot  produce 

t 

a  continuous  minority  carrier  current.  For  this  reason,  an  additional 
component  of  minority  carrier  current  arises  through  the  mechanism  of 
thermal  diffusion;  in  combination,  drift  and  diffusion  produce  electric  , 
current  continuity  (due  to  minority  carriers)  within  a  reverse  biased 
high-low  semiconductor  junction. 

Figures  6  and  7  show  the  calculated  majority  carrier  and  minority 
carrier  distributions,  respectively,  in  a  reyerse  biased  hiph-low 
junction.  This  sequence  of  calculations  demonstrate  that  minority 
carrier  depletion  takes  place  within  both  sides  of  this  structure,  and  1 
that  the  amount  of  depletion  increases  with  an  increase  of  electric 
current.  Throughout  these  regions  of  carrier  depletion  charge  neutrality 
is  maintained  by  the  removal  of  an  almost  equal  number  of  minority  and 
majority  carriers,  although  the  change  of  majority  carriers,  due  to  this 
mechanism,  is  too  small  to  be  observed  in  Fig.  6.  .  1 

i 

1.3  Calculated  JFET  Characteristics  i 

Figure  8  presents  a  comparison  between  the  calculated  volt-ampere 
characteristics  of  the  two  JFET  structures  shown  in  Figs.  1  and  Z\  these 
calculations  are  in  general  agreement  with  tM>$e  presented  by  M.  Reiser. 

In  particular,  it  is  shown  here  that  a  substantia!  decrease  of  saturation 
resistance  (assuming  identical  channel  dimensions  and  channel  doping)  can 
be  expected  in  a  JFET  containing  a  low  conductivity  semiconductor  substrate 
and  a  field  dependent  carrier  mobility  within  its  source-drain  channel. 

In  contrast,  Fig.  8  shows  this  difference  of  ‘aturation  resistance  is 
significantly  reduced  when  a  constant  carrier  mobility  is  assumed  in 
these  structures.  From  Fig.  8,  it  is  clear  that  this  reduction  of 

.  ’  i 

saturation  resistance  is  attributable  to  the  low  conductivity  semiconduc¬ 
tor  substrate,  in  conjunction  with  the  assumption  of  a  field  dependent 
carrier  mobility. 

Figures  9  and  10  illustrate  the  calculated  mobile  carrier  distribution 
within  this  JFET  containing  a  substrate  (Fig.  2),  when  the  structure  is 
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.  Calculated  minority  carrier  distribution  in  a 
reverse  biased  high-low  junction  (silicon). 
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Fig.  8.  Calculated  volt-ampere  charcteristics  cf  the 
JFET  structures  shown  in  Figs.  1  and  2. 


Fig.  10.  Calculated  contours  of  c' w*nt  majority  carrier  density  within 
a  JFET  containing  a  su  .  (field  dependent  carrier  mobility). 


biased  well  Into  current  saturation.  Shown  In  this  illustration  are  < 
contours  of  constant  electron  density  as  measured  with  respect  to  the 
Impurity  atom  density  in  the  source-drain  channel  (1.5  x  10;  atoms^cm  ); 

for  example,  the  10X  contour  represents  a  mobile  electron  of  1.5  x 

16  3  1  1 

10  electrons/cm  .  Figure  9  is  based  upon  an  assumption  that  mobile 

carriers  exhibit  a  constant  mobility  throughout  the  entire  structure. 

Figure  10  shows  the  calculated  mobile  electron  distribution  when  these. 

carriers  are  assumed  to  exhibit  a  field  dependent  drift  mobility.  1 

From  Fig.  9,  a  constant  carrier  mobility  within  the  source-drain 

•  j 

channel  produces  a  situation  where  the  substrate  carrier  distribution 

i 

differs  little  from  the  distribution  implied  by  a  one-dimensional 
analysis  of  high-low  junctions  (Fig.  3).  A  layer  of  excess  »wbile  car- 

i 

riers  lies  essentially  parallel  to  the  metallurgical  high-low  junction, 
and  this  layer  exhibits  a  small  contraction  near  the  drain  side  of  the  1 
structure.  A  detailed  calculation  of  the  electric  current  distribution 
near  the  "pinched-off"  region  shows  that  a  substantial  part  of  the1  total 

i 

source-drain  current  arises  from  carrier  transport  thpough  this  layer  of 
carrier  accumulation  For  all  practical  purposes,  the  effective  channel 
width  is  increased  by  this  high  conductivity  layer  residing  at  the 
channel -to-substrate  boundary f 

2  4  ' 

In  a  previous  investigation  *  it  was  shown  that  electric  currenc 

I  i 

continuity  is  maintained  in  a  JFET  (Fig.  1)  by  the  process  of  carrier 
accumulation. '  If  the  total  electric  current  entering  the  "pinched-off" 
region  exceeds  the  total  electric  current  within  this  "pinched-off"  region 
(due  to  velocity  limited  carrier  transport)  carrier  accumulation  will  take 
place,  and  this  accumulation  provides  the  means  whereby  electric  current 
continuity  Is  maintained  between  the  source  and  drain.  Furthermore,  it 
was  shown  that  this  region  of -carrier  accumulation  will  physically  expand 
the  conduction  path  through  this  "pinched-off"  region,  and  even  produce 
contraction  of  the  gate  junction  space-charge  layer.  It' is  suggested 
here  that  the  calculated  substrate  carrier  accumulation  shown  in  Fig.  10 
(and  also  in  the  publication  of  M.  Reiser^)  is  another  manifestation  of 
accumulation  region  expansion. 


For  this  reason,  the  soft  saturation  characteristics  attributed 
to  the  substrate  in  a  JFET  (Fig.  2)  appears  to  result  from  an  additional 
degree  of  freedom  for  carrier  accumulation  within  the  source-drain  channel 
(when,  compared  with  the  structure  ^hown  in  Fig.  1);  In  the  absence  of  a 
substrate,  carrier  accumulation  is  the  dominant  mechanism  whereby  current 
continuity  Is  maintained;  expansion  of  this  region  of  accumulation  Is 
limited  by  tne  gate  junction  space-charge  layer  and  the  Insulating  sur- 

i  '  :  ; 

face.  In  contrast,  structures  containing  a  substrate  provide  semiconduc¬ 
tor  material  in, which  this  region  of  excess  carriers  can  expand  and 

thereby  .produce  a  soft  current  saturatibn  characteristic. 

1  . 

) 
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CHAPTER  II 


The  Insulated  Gate  Field-Effect  Transistor 
2.0  Introduction 

Early  Investigations  of  IGFET  operation^  suggest  that  electric  cur¬ 
rent  saturation  Is  attributable  to  source-drain  channel  pinch-off  in  a 

manner  similar  to  that  proposed  by  Shockley  in  his  publication  on  the 

2 

unipolar  transistor.  Briefly,  it  was  suggested  that  current  satura¬ 
tion  Is  a  consequence  of  the  potential  distribution  between  two  regions 
within  the  source  drain  channel: 

Region  I  —  from  the  source  junction  to  the  channel  pinch-off 
point,  and 

Region  II  —  from  the  channel  pinch-off  point  to  the  drain 
junctl on. 

This  concept  of  IGFET  operation  has  been  elaborated  on  by  many  authors, 
yet  there  appears  to  be  little  understanding  concerning  the  detailed 
mechanisms  contributing  to  this  potential  distribution. 

Another  viewpoint  concerning  electric  current  saturation  is  that 
channel  pinch-off  produces  a  region  of  small  electrical  conductivity 
within  the  source-drain  channel,  and  this  region  of  small  conductivity 
produces  the  observed  saturation  resistance.  Detailed  mechanisms 
involved  in  channel  pinch-off  have  been  studied  by  many  workers.  For 

4 

example,  Chiu  and  Sah  present  a  two-dimensional  analytical  solution 
for  the  potential  distribution  throughout  this  region;  due  to  inherent 
limitations  of  their  analysis  it  provides  little  insight  into  the 

5 

solution  of  this  problem.  Further,  Pao  and  Sah  view  channel  pinch-off 
as  a  "bottleneck"  in  which  diffusion  becomes  an  important  mechanism  of 
carrier  transport.  Two-dimensional  numerical  solutions  of  IGFET  operation 
were  presented  by  M.  B.  Barron®  and  by  Vandorpe  and  Xuong;^  both  of  these 
numerical  solutions  show  the  <  xi  cnc'.  r  o  rl!  defined  channel 
pinch-off  point,  althouth  it  has  been  established  that  this  result  is 
a  consequence  of  computational  errors  associated  with  their  numerical 

g 

method  (see  Chapter  IV). 

The  purpose  of  this  chapter  is  to  present  the  initial  results  of 
a  detailed  two-dimensional  steady-sto.^  mathematical  investigation  of 
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IGFET  operation.  The  system  of  equations  applied  to  this  problem  were 

g 

previously  outlined  by  Van  Roosbroeck  in  connection  with  his  studies 
of  the  transport  of  electrons  and  holes  in  a  semiconductor.  Finite 
difference  methods  are  used,  with  the  aid  of  an  electronic  computer. 

Using  this  technique,  a  rigorous  analysis  is  under  way  with  an  aim 
toward  understanding  the  physical  mechanisms  of  operation  that  produce 
the  electrical  characteristics  exhibited  by  an  insulated-gate 
field-effect  transistor. 

This  mathematical  technique  introduces  many  difficulties  of  inter¬ 
pretation.  The  system  of  equations  used  In  our  two-dimensional  computer 
solutions  of  the  IGFET  problem  are  exceedingly  general,  and  implicitly 
contain  all  mechanisms  Important  to  the  operation  of  this  semiconductor 
device.  As  a  consequence,  it  is  sometimes  difficult  to  explain  the 
results  of  these  computer  calculations  in  terms  of  recognized  physical 
mechanisms  involved  in  the  operation  of  an  IGFET.  One  method  for 
solving  this  difficulty  Is  to  quantitatively  and  qualitatively  compare 
the  results  of  our  computer  calculations  with  elementary  theory  taken 
from  textbooks  and  technical  papers;  this  elementary  theory  Is  based 
upon  specific  physical  mechanisms  of  device  operation.  If  the  computer 
calculations  are  in  agreement  with  elementary  theory  it  is  presumed  the 
dominant  physical  mechanisms  are  known.  If,  instead,  computer  calculations 
disagree  with  elementary  theory,  it  is  presumed  that  the  elementary  theory 
is  Incomplete;  seldom  can  it  be  shown  that  elementary  theory  is  incorrect. 

After  establishing  that  a  particular  aspect  of  the  existing  elementary 
theory  is  inadequate,  an  important  technique  for  determining  the  mechanisms 
involved  in  semiconductor  device  operation  (and  a  technique  used  exten¬ 
sively  in  our  investigations)  is  to  conduct  experiments  on  the  computer. 
These  experiments  may  consist  of  new  or  unusual  geometrical  configurations, 
biasing  conditions,  etc.;  each  experiment  designed  to  simplify  the 
problem  of  interpretation.  In  this  fashion,  we  can  often  determine  the 
individual  mechanisms  involved  in  IGFET  operation  and,  hopefully,  put 
together  an  elementary  theory  that  is  in  substantial  agreement  with  the 
rigorous  computer  calculations.  The  overall  purpose  of  this  effort  Is 
to  evaluate  existing  theoretical  concepts  of  IGFET  operation  and,  where 


necessary,  modify  and  improve  these  concepts  to  yield  an  elementary 
(or  engineering)  theory  that  can  be  used  without  a  computer. 

2.1  Analysis 

It  has  been  shown  that  the  hole  and  electron  distributions  in 

q 

semiconductor  material  are  described  by  the  equations 

a)  div  grad  ^  ^  (N  -  n  +  p) 

KEq 

b)  Jp  =  -  qDp  grad  p  -  qupp  grad  * 

c)  Jn  a  qDn  grad  n  -  qnnn  grad  (1) 


d)  div  Jp  =  qRp 

e)  div  Jn  =  q^  and 

f)  =  ^p  +  ^n 

when  it  is  assumed  that  no  trapping  mechanisms  exisc  within  the 
structures  under  consideration.  Finite  difference  methods  are  used 
to  numerically  solve  this  system  of  equations;  details  of  these 
computational  methods  are  outlined  in  Chapter  IV. 

As  in  all  mathematical  investigations  of  this  type  there  are 
important  physical  parameters  associated  with  IGFET  operation  for 
which  little  (or  no)  information  is  available.  For  example,  little 
quantitative  date  is  available  concerning  either  the  small  field  or 
large  field  mobility  of  carriers  in  a  source-drain  channel.  Similarly, 
minority  carrier  lifetime  is  seldom  known  with  any  degree  of  accuracy. 

A  consequence  of  this  situation  is  the  necessity  to  make  assumptions 
concerning  the  physical  parameters  of  an  IGFET  that  i.iay  not  be  com¬ 
pletely  correct.  Throughout  this  mathematical  study  we  have  adopted 
values  for  specific  physical  p«ra ..^ters  that  are  consistent  with 
information  presently  available  in  the  technical  literature.  In 
most  situations  the  quantitative  values  adopted  for  these  parameters 
have  little  infuence  upon  the  physical  mechanisms  under  investigation; 
when  this  influence  is  significant,  the  topic  is  discussed  in  detail. 


It  has  been  suggested  that  the  large  field  transport  charac¬ 
teristics  of  mobile  carriers  within  a  source-drain  channel  exhibit 
hot  electron  characteristics  that  differ  from  those  encountered  in 
bulk  semiconductor  material.10  This  suggestion  appears  to  be  con- 

ii  13 

firmed  by  recent  publications  on  this  topic.  From  these 

publications ,  we  have  assumed  in  the  present  analysis  that  mobile 

carriers  exhibit  a  maximum  (or  terminal)  velocity  of  about 

8.6  x  10^cm/sec;  this  magnitude  is  taken  from  measurements  upon 
14 

bulk  silicon,  and  it  is  within  a  factor  of  2  of  the  measured 

13 

terminal  velocity  for  carriers  within  an  MOS  channel.  From  these 
publications  we  have  alsc  used  traditional  assumptions  for  the  warm 
electron  region  of  the  carrier  velocity:12*1^ 

v  a  E  when  0  <  v  <  v„ 

—  c 

1/2 

v  a  E  '  when  v„  <  v  <  vm  , 
c  m 

where  E  represents  the  electric  field  and  the  symbols  v,  vc-  and  vm 
represent  the  carrier  velocity,  the  critical  carrier  velocity  producing 
a  relation  v  a  E1/2,  and  the  maximum  (or  terminal)  carrier  velocity, 
respectively. 

The  magnitude  of  electric  field  producing  this  critical  carrier 

velocity  (vj  remains  an  unanswered  question.  It  was  recently 

llc12  3 

suggested  *  that  an  electric  field  of  1.5  x  10  cm/sec  produces 

warm  electrons  (v  a  E1^2)  in  an  MOS  channel  exhibiting  a  low  field 

2 

channel  mobility  of  about  250  cm  /volt  sec.  Therefore,  we  have 
adopted  this  value  of  critical  carrier  velocity  in  the  present 
two-dimensional  analysis,  and  the  calculated  results  based  upon 
this  assumption  appear  to  be  in  agreement  with  qualitative  measure¬ 
ments  upon  laboratory  MOSFET  structures. 

Our  present  method  for  solving  Eqs.  1  is  based  upon  an  assumption 
that  no  carrier  recombination  exists  within  the  structure  under  con¬ 
sideration.  It  is  believed  the  small  recombination/generation  rate 
present  within  modern  silicon  IGFET  structures  has  little  influence 
upon  its  steady-state  electrical  characteristics.  Because  the  present 
investigation  is  restricted  to  the  steady-state  properties  of  this 
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semiconductor  device,  there  is  no  reason  to  suspect  that  by  neglecting 
recombination/generation  in  our  mathematical  model  we  have  introduced 
a  significant  computational  error. 


Throughout  this  mathematical  investigation,  particular  emphasis 
is  placed  upon  the  selection  of  boundary  conditions  that  do  not  intro¬ 
duce  errors  in  the  calculated  results.  In  particular,  the  boundary 
conditions  used  in  this  analysis  approximate  physical  and  electrical 
characteristics  at  the  outer  periphery  of  this  semiconductor  structure, 
rather  than  along  internal  boundaries  established  by  its  physical  or 
electrical  properties.  A  homogeneous  impurity  atom  distribution  is 
assumed  within  both  the  n-type  and  p-type  semiconductor  material.  The 
exposed  semiconductor  surface  is  assumed  to  be  an  ideal  electrical 
insulator;  i.e.,  no  electric  current  is  permitted  normal  to  these 
bounding  surfaces.  The  ohmic  contacts  are  approximated  by  equipotential 
surfaces  that  are  charge  neut.ral ;  these  ohmi  c  contacts  have  been 
located  sufficiently  far  from  the  active  regions  of  this  device  to 
have  no  Influence  upon  its  electrical  properties. 


In  the  calculations  presented  here,  the  physical  and  geometrical 
properties  of  our  analytical  model  (Fig.  11)  were  selected  as  represen¬ 
tative  of  a  typical  MOS  field-effect  transistor.  All  calculations  are 
based  upon  analytical  models  that  have  one  of  three  different  distances 
between  the  source  and  drain  metallurgical  junctions;  10  Mm,  5.0  pm, 
and  1.0  pm.  Rather  than  calculate  the  source  and  gate  junction  shapes 
shown  in  Fig.  11,  the  simplified  model  shown  in  Fig.  12  is  used  to 
approximate  an  IGFET;  geometrical  differences  between  Figs.  11  and  12 
have  little  Influence  upon  the  conclusions  presented  here. 


In  the  source  and  drain  regions  of  these  structures  a  donor  atom 
19  3 

density  of  10  atoms/cm  is  assumed  to  form  ideally  abrupt  p-n  junctions 
into  a  substrate  region  containing  ai.  acceptor  atom  density  appropriate 
to  the  channel  length; 

Channel  length  (pm)  Substrate  dop:ng  ( atoms /cm  ) 

10  2  x  1015 

5.0  2  x  1015 

1.0  2  x  1016  . 
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Fig.  12.  Computational  model  used  In  this  1nvest1;at1on 


To  obtain  a  reasonable  comparison  between  the  electripal  properties 
of  these  fcitalytical  models,  it  was  necessary  to  introduce  modifications 
to  compensate  for  changes  of  substrate  doping,  and  thereby  obtain  a! 
consistent  threshold  voltage.  For  this  reason,  an  oxide  (SiO«)  thick- 

O  ^ 

ness  of  1000  A  is  assumed  in  conjunction  with  the  substrate  doping  of 
2  x  10^5atoms/cm3  and  288  A  with  the  doping  of  2  x  10^6atoms/cm3.  ! 

2.2  Source-Drain  Channel  “Pinch-off11  in  an  IGFET  i  ,  i 

i 

2 

In  his  investigations  of  JFET  operation,. Shockley  showed  that 

channel  pinch-off  by  the  gate  junctions  produces  a  region  of  carrier 

depletion  across  which  appears  that  nortion  of  the  source-drain  voltage 

in  excess  of  the  channel  pinch-off  voltage.  Because  current  continuity 

must  be  maintained  throughout  this  semiconductor  structure,  it  is  clear 

that  channel  pinch-off  in  a  JFET  cannot  produce  complete  depletion  of 
15 

mobile  carriers;  a  sufficient  carrier  density  will  always  exist  , 
throughout  the  pinch-off  region  to  maintain  the  source-drain  electric 
current.^6  A  similar  statement  can  be  made  for  the  pinch-off  region 
of  an  IGFET  if,  indeed,  the  mechanisms  of  operation  in  an  IGFET  are  ' 

equivalent  to  those  '■"■countered  in  a  JFET..  The  purpose  of  this  section 
is  to  consider  the  .aical  mechanisms  encountered  in  thei source-drain 
channel  of  an  IGFET  v  n  it  is  assumed  that  the  structure  is  biased  well 
into  electric  current  i.  ^ration. 

I 

For  purposes  of  simplicity,  it  is  assumed  here  that  the  source- 

i 

drain  channel  pinch-off  point  in  an  IGFET  is  coincident  with  that 
location  where  the  difference  of  potential  between  the  semiccnductor- 
insulator  interface  and  the  gate  electrode  just  equals  the  threshold 
voltage.  A  short  distance  from  this  location  the  gate  electrode  and 
the  semiconductor  are  at  the  same  potential  and  no  electric  field 
exists  in  a  direction  perpendicular  to  the  semiconductor' surface, 

Fig.  13.  On  the  source  side  of  this  location  the  gate  induced  electric  ■ 
field  produces  an  accumulation  of  minority  carriers  along  the  semicon¬ 
ductor  surface.  In  contrast,  near  the  drain  side  of  this  location  the  gate 
induced  electric  field  forces  mobile  carriers  away  from  the  semiconductor 
surface;  therefore,  between  the  channel  pinch-off  point  (approximately)  and 
the  drain  the  semiconductor  surface  is  essentially  free  of  mobile  carriers. 

i 

This  situation  has  a  significant  influence  upon  the  mechanism  of  chan¬ 
nel  pinch-off.  Between  the  drain  junction  and  the  channel  pinch-off  point 


the  difference  of  potential  is  established  almost  entirely  by  ionized  im¬ 
purities  within  the  drain  junction  space-charge  layer.  Further  because 
this  difference  of  potential  is  always  smaller  than  the  applied  drain- 
source  voltage,  the  pinch-off  point  must  reside  inside  this  drain  junction 
space-charge  layer,  a  region  of  large  longitudinal  electric  field. 

If,  indeed,  channel  pinch-off  produces  a  "bottle-neck"  mechanism 
for  the  source-drain  electric  current  (similar  to  that  encountered  in 
a  JFET)  a  constraint  must  exist  upon  the  physical  width  of  this  con¬ 
ducting  channel.  In  a  JFET  this  constraint  is  a  consequence  of  an 
electric  field  arising  from  the  gate  junction  space-charge  layers.^6  In 
contrast,  the  channel  pinch-off  point  in  an  IGFET  is  located  within  a 
region  of  large  longitudinal  electric  field;  this  field  provides  a 
mechanism  whereby  carriers  can  (and  do)  leave  the  channel  and  continue 
toward  the  drain  junction.  For  this  reason,  pinch-off  of  the  source- 
drain  channel  does  not  produce  a  constriction  of  conducting  carriers 
at  the  semiconductor  surface  but,  instead,  it  produces  an  increase  in 
the  width  of  this  conducting  channel.  It  is  therefore  suggested  that 
pinch-off  of  mobile  carriers  will  not  exist  in  a  structure  biased  well 
into  electric  current  saturation.  This  qualitative  argument  has  been 
verified  by  a  rigorous  numerica7  calculation  of  the  mobile  carrier  and 
electric  flux  distributions  in  an  IGFET. 

To  illustrate  this  situation,  Figs.  14  and  15  show  the  calculated 
minority  carrier  distributions  in  an  IGFET  biased  below  saturation  and 
in  one  biased  well  into  current  saturation,  respectively.  Below  sat¬ 
uration  (Fig.  14)  minority  carriers  form  a  definite  channel  between  the 
source  and  drain.  In  contrast,  when  this  structure  is  biased  into  cur¬ 
rent  saturation  (Fig.  15)  mobile  carriers  start  to  move  away  from  the 
semiconductor  surface  at  a  location  far  removed  from  what  is  normally 
considered  the  pinch-off  point. 

From  Fig.  15,  it  is  difficult  to  maintain  the  concept  of  channel 
pinch-off  in  an  IGFET.  This  situation  is  further  emphasized  in  the  cal¬ 
culated  mobile  carrier  flux  distribution  shown  in  Figs.  16  and  17.  Figure 
16  sho^s  the  flux  distribution  in  an  IGFET  biased  well  below  current  satur¬ 
ation.  Similarly,  Fig.  17  shows  the  flux  distribution  when  this  structure 
is  biased  well  into  current  saturation.  Clearly,  little  is  found  to  justify 
the  existence  of  pinch-off  condition  in  the  theory  of  IGFET  operation  (see 
Fig.  17).  Therefore,  throughout  this  report  the  traditional  concept  of  a 
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Fig.  14.  Calculated  minority  carrier  distribution  in  an  MOSFET 
assuming  a  constant  carrier  mobility  (V^g*  3.0  volts j 
Vcn3  1 .0  volt) . 
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Fig.  15.  Calculated  minority  carrier  distribution  in  an  MOSFET 
assuming  a  constant  carrier  mobility  (VSG=  3.0  volts; 
V<;.n=  6.0  volts). 


Fig.  16.  Calculated  mobile  carrier  flux  distribution  in  an  MOSFET 
assuming  a  constant  carrier  mobility  (Vcr=  3.0  volts; 

Vq n=  1.0  volt).  SG 


Oxidt 


I 


I 


31 


i 


E 


'  '  !  J  i 

I 

1  5  f 

)  ’  I  i 

channel  pinch-off  point  is  redesignated  as  a  point  of,  channel  termination; 
thereby  we  are  eliminating  any  implication .that  channel  pinch-off  produces1 
a  region  of, large  electrical  resistance,  as  encountered  in  a  OFET. 

■  •  1  1  I  ,  , 

2.3  Electric  Currerit  Saturation 

Fjrom  our  calculations  ,of  IGFET  operation  itiis  concluded  that 
two  distinctly  different  mechanisms  are  capable  of  producing  electric 

!  current  saturation.  In  a  very  short  channel  device  electric  current 
saturation  can  result  from  velocity  saturation  of  mobile  carriers  within 
the  source-drain  channel.  In  a  long-channel , IGFET  velocity  saturation 

i  '  i  |  '  1  1 

{  of  mobile  carriers  may  take  plac4,  although  this  mechanism  has  little 
influence  upon  its  volt-ampere  characteristics.  ,Electric  current  sat¬ 
uration  in  a  long-channel  IGFE;T  is  a  consequence  of  the  potential 
distribution  associated  with  the  drain  junction  space-charge  layer. 

i  *  1  ,  i 

1  It  is  interesting  to  note,, both  experiment  and  theory  show  that 
electric  current  saturation  can  take  place  in  an  IGFET  containing 'no 

1  i  j  1 

velocity  saturation,  and  when  operating  under  circumstances  that  do  not 
permit  channel  "pinch-off",.  \ 

I 

In  a  long-channel  IGFET  electric  current  saturation  takes  place 
when  the  applied  source-drain  voltage  is  sufficient  to  deplete  mobile 
carriers  from  the  drain  end  of  the  channel.  One-dimensional  theory 
shows  this  depletion  can  arise  in  two  different  ways:  first,  by  a 
1  reduction  of  voltage  between  the  gate  electrode  and  the  semiconductor 

surface  and,  second,1  by  an  increase  of  voltage  between  the  semiconductor 
surface  and  charge  neutral  regions  of  the  substrate.  A  reduction  of 

1  gate  to  surface  voltage  produces  a  reduction  of  electric  field  normal 
to  the  semiconductor  surface;  thereby  a  reduction  is  obtained  in  the 
total  surface  charge  within  the  semiconductor  (Q^  +  Q0).  Assuming  the 
gate  voltage  is  above  threshold  for  the  device  under  consideration 

«  I 

VSg),  this  reduction  of  electric  field  produces  an  equivalent 

reduction  of  mobile  charges  in  the  inversiori  layer  (Q..).  Similarly, 

one-dimensional  theory  also  shows  that  if  a  given  electric  field 

exists  at  the  semiconductor  surface  (hence,  Q.  +  QD  is  specified) 

an  increase  of  surface  potential  (as  measured  with  respect  to  charge 

neutral  regions  of  the  substrate)  produces  an  increase  of  depiction 

charge  (Qp)  and,  therefore,  a  reduction  of  mobile  charges  in  the  , 

inversion  layer  (Q. ). 

32  1 


i 


1 


Deep  saturation  takes  place  In  a  long-channel  IQFET  when  the 
source-drain  voltage  exceeds  the  source-gate  voltage.  This  mode  of 
operation  produces  a  situation  where,  at  some  location  near  the  drain 
junction,  the  gate  induced  electric  field  becomes  zero  and,  hence,  the 
source-drain  channel  is  terminated  (see  Figs.  13  and  17).  Between 
this  channel  termination  point  and  the  drain  a  negligible  carrier 
density  exists  along  the  semiconductor  surface;  throughout  this 
region  the  longitudinal  potential  distribution  is  determined  by 
uncompensated  ionized  impurities  associated  with  the  drain  junction 
space-charge  layer.  Therefore,  any  change  of  drain  voltage  (assuming 
it  remains  larger  than  the  gate  voltage)  produces  a  change  of  voltage 
between  the  point  of  channel  termination  and  the  drain,  and  the 
voltage  applied  to  the  source-drain  channel  remains  essentially  con¬ 
stant.  This  mechanism  was  first  postulated  by  IhantolaJ  and  has 
been  verified  by  rigorous  computer  calculations  of  this  semiconductor 
device. 

Figure  18  shows  the  calculated  volt-ampere  characteristics  of  a 
typical  long-channel  IGFET.  In  this  particular  calculation  It  is  assumed 
that  mobile  carriers  within  the  source-drain  channel  exhibit  no  velocity 
saturation  (constant  drift  mobility).  From  Fig.  18,  it  is  evident  that 
current  saturation  is  initiated  at  a  source-drain  \oltage  that  is  below 
the  value  producing  channel  termination.  Tnis  initiation  of  current 
saturation  (at  a  reduced  drain  voltage)  is  assumed  to  result  from  the 
onset  of  channel  termination,  although  from  previous  discussions  it  is 
not  obvious  how  this  situation  arises  in  an  IGFET  structure. 

From  one-dimensional  IGFET  theory,  the  total  electrostatic  charge 
at  the  semiconductor  surface  (Qy)  Is  attributable  to  uncompensated 
ionized  impurity  atoms  within  a  region  of  majority  carrier  depletion 
(Qd)  and  mobile  carriers  within  a  t'r>  inversion  layer  (Q^).  This 
total  charge  (Qy)  must  be  proportional  to  the  normal  component  of 
electric  field  at  the  semiconductor  surface, 
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Fig.  18.  Calculated  volt-ampere  characteristics  of  a 
long-channel  IGFET  assuming  a  constant  drift 
mobility  for  carriers  in  the  source-drain 
channel.  (1^=1000  ft,  L  =10  ym). 


assuming  the  structure  is  free  of  interface  charges.  In  addition,  the 
total  depletion  charge  (Q^)  is  assumed  to  be  determined  by  the  surface 
to  substrate  voltage,  giving 


Qg  =  ^dxD  » 

where  xQ  (the  depletion  layer  width)  is 

o  2kc£aV. 

Y2  -  s  0  S 

*D  • 

Hence,  the  inversion  layer  charge  (Q,.)  is  proportional  to 

Qi- 


(3) 

(4) 

(5) 


From  eq.  5,  it  is  implied  that  a  variation  of  surface  voltage 
(Vs),  holding  (VQ  -  V$)  constant,  will  produce  a  variation  of 
inversion  layer  charge  (Q^).  In  two-dimensional  IGFET  theory  this 
conclusion  is  not  obvious.  Poisson's  equation  in  two  dimensions 
shows  that  a  change  of  voltage  between  the  semiconductor  surface  and 
neutral  substrate  material  does  not  necessarily  imply  a  change  of 
depletion  layer  charge;  a  redistribution  of  charge  could  arise  to 
support  this  change  of  voltage.  Specifically,  in  some  regions 
ionized  impurities  are  shared  between  the  gate  electrode  and  the 
drain  (i.e.  some  are  electrostatically  tied  to  the  gate  and  others 
to  ionized  impurities  in  the  drain  region).  In  two  spatial 
dimensions,  only  those  impurities  electrostatically  tied  to  the 
gate  electrode  contribute  to  the  surface  potential.  For  this 
reason,  a  redistribution  of  electrostatic  association  between  ionized 
impurities  in  the  depletion  layer,  the  gate  electrode,  and  the  drain 
contact  could  arise  with  no  change  of  the  depletion  charge  Qg  in 
eq.  5. 
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In  order  to  evaluate  this  situation  the  volt-ampere  characteristics 

have  been  calculated  for  an  IGFET  containing  an  unusually  large  oxide 

thickness  (Fig.  19).  This  large  oxide  thickness  makes  it  necessary  in 

these  calculations  to  assume  a  source-gate  bias  in  excess  of  about 

100  volts  before  a  significant  conducting  channel  is  formed  at  the 

semiconductor  surface.  From  Fig.  20,  at  a  source-gate  bias  of  108  volts 

this  semiconductor  device  exhibits  electric  current  saturation  when  the 

source-drain  bias  is  in  excess  of  about  2.0  volts;  therefore  (Vq-Vr) 

in  eq.  5  is  nearly  constant.  Although  there  is  little  (or  no) 

experimental  information  on  an  IGFET  of  this  type,  other  situations 

arise  in  the  operation  of  more  conventional  devices  that,  to  all 

17 

practical  purposes,  experimentally  vefify  our  calculations. 

The  calculated  electrical  properties  of  this  wide  gate  IGFET 
(Fig.  20)  exhibit  electric  current  saturation,  despite  the  fact  that 
an  essentially  constant  gate  induced  electric  field  exists  along  the 
semiconductor  and  insulator  interface.  From  traditional  theory  of 
IGFET  operation,  current  saturation  in  this  structure  is  attributable 
to  an  increase  of  surface  potential  arising  from  the  application  of  a 
drain  voltage.  An  increase  of  surface  voltage  produces  an  increase  of 
depletion  charge  (Q0)  and,  hence,  a  decrease  of  Inversion  charge  (Q^); 
this  decrease  of  inversion  charge  along  the  insulator  and  semiconductor 
interface  is  shown  in  Figs.  21  and  22.  A  rigorous  calculation  of  this 
inversion  layer  carrier  distribution  is  qualitatively  consistent  with 
one-dimensional  theory,  although  it  is  necessary  to  verify  that  a 
decrease  of  is  equal  in  magnitude  to  an  increase  of  depletion  layer 
charge  (QD). 

From  a  detailed  study  of  the  calculated  electric  field  distribution 
near  the  inversion  layer  of  this  device,  it  has  been  qualitatively 
established  that  the  inversion  charge  (Q^)  and  the  depletion  charge  (QD) 
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le  carrier  distribution  in  the  source-drain 
GFET  containing  a  wide  insulator  (Fig.  19). 


39 


satisfy  the  equality 


QT  =  Ql  +  QD  »  (6) 

where  Is  the  total  charge  at  the  semiconductor  surface.  This  is  not 

Intended  to  imply  that  our  calculations  are  in  agreement  with  traditional 

1  ’ft 

concepts  of  IGFET  operation.  For  example,  traditional  theory  *  shows 
that  this  wide  insulator  structure  will  exhibit  channel  termination  (or 
pinch-off)  when  the  drain  voltage  is  about  3.2  volts  and  the  gate  voltage 
is  about  154  volts.  Figure  22  shows  substantial  differences  between  our 
computer  calculated  results  and  the  conclusions  drawn  from  this  tra¬ 
ditional  theory. 

From  Flo.  22,  the  minimum  carrier  density  (V^  =  4)  is  about  80 
times  the  background  doping;  this  magnitude  is  substantially  larger 
than  one  would  estimate  from  present  IGFET  theory.  Furthermore,  at  a 
drain  bias  of  6.0  volts  the  point  of  minimum  inversion  carrier  density 
remains  close  to  the  drain  junction  (Fig.  22),  rather  than  at  a 
location  where  the  surface  voltage  is  approximately  4.0  volts.  This 
situation  shows  that  some  fundamental  discrepancies  exist  between  our 
traditional  theory  of  IGFET  operation  and  the  results  obtained  from  a 
rigorous  solution  of  the  problem. 

Substantial  effort  has  been  directed  toward  the  development  of 
more  accurate  approximations  for  the  mechanisms  involved  in  electrical 
conduction  through  an  inversion  layer.  There  are  many  obvious 
simplifications  in  the  available  one -dimensional  theory  (gradual  channel 
approximation,  neglecting  diffusion,  etc.)  although,  thus  far,  by 
eliminating  these  simplifications  we  have  only  complicated  the  equations, 
with  no  fundamental  improvement  in  the  results.  From  this  effort',  it  is 
concluded  that  important  mechanisms  are  missing  from  our  present  one- 
dimensional  theory  of  electrical  conduction  in  a  source-drein  channel, 
yet  we  have  thus  far  been  unable  to  identify  these  mechanisms. 

Calculations  show  that  velocity  saturation  in  a  long-channel  IGFET 
is  restricted  to  the  drain  junction  space-charge  layer.  Through  a  com¬ 
parison  of  Figs.  23  and  24,  it  is  concluded  that  velocity  saturation  in 
this  region  of  the  device  has  a  negligible  influence  upon  its  volt-ampere 
characteristics.  The  reasons  for  this  situation  are  evident  from  the  • 
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calculated  mobile  carrier  distribution,  particularly  in  the  immediate 
vicinity  of  the  drain.  Briefly,  velocity  saturation  can  do  little  to 
influence  the  electric  current  if,  for  all  practical  purposes,  the 
total  number  of  carriers  producing  this  current  can  increase  indefinitely. 
In  this  situation,  a  reduction  of  drift  mobility  (with  an  increase  of 
electric  field)  results  in  an  increase  of  carriers  and,  thereby,  current 
continuity  is  maintained  between  the  source  and  drain. 

Figures  25  and  26  illustrate  this  mechanism,  as  calculated  for  a 
long-channel  IGFET.  A  comparison  between  Fig.  15  (constant  drift 
mobility)  and  Fig.  25  (field  dependent  drift  mobility)  shows  that  by 
•introducing  a  field  dependent  mobility  into  these  calculations  a  sub¬ 
stantial  increase  is  observed  in  the  number  of  mobile  carriers  residing 
within  the  drain  junction  space-charge  layer.  Further,  if  for  purposes 
of  experiment  a  decrease  is  introduced  into  the  maximum  carrier  velocity, 
a  proportional  increase  is  observed  in  the  number  of  carriers  stored 
within  this  region  of  the  structure.  As  previously  stated,  this  increase 
of  carrier  density  in  proportion  to  a  decrease  of  carrier  mobility  is 
interpreted  to  represent  the  mechanism  whereby  velocity  saturation  has 
a  negligible  influence  upon  the  volt-ampere  characteristics  of  an  IGFET. 

This  mechanism  of  spreading  within  the  drain  space-charge  region 
is  also  observed  in  the  calculated  electron  flux  distribution.  Fig.  26. 

A  comparison  between  Fig.  17  (constant  drift  mobility)  and  Fig.  26 
(field  dependent  drift  mobility)  shows  a  substantial  spread  of  electron 
flux  throughout  the  drain  space-charge  layer,  with  the  introduction  of 
a  field  dependent  drift  mobility. 

In  contrast  with  the  long-channel  structure  (Lc  =  10  um), 

velocity  saturation  in  a  short-channel  IGFET  (L  =1.0  ym)  has  a 

c 

significant  influence  upon  its  volt-ampere  characteristics.  Figure  27 
shows  that  in  the  absence  of  velocity  saturation  a  short  channel  IGFET 
will  exhibit  a  relatively  soft  electric  current  saturation.  A  decrease 
of  channel  length,  with  an  increase  of  drain  voltage,  is  traditionally 
considered  the  mechanism  producing  soft  current  saturation,  as  observed 
in  short-channel  devices,  although  this  mechanism  is  probably  of  little 
consequence.  In  an  IGFET  of  these  dimensions  (Lc  =  1.0  ym),  velocity 
saturation  produces  a  sufficient  modification  of  the  volt-ampere  charac¬ 
teristics  to  probably  mask  any  consequences  of  changes  in  channel  length. 


Figure  28  demonstrates  the  Influence  of  velocity  saturation 

upon  the  volt-ampere  characteristics  of  this  short-channel  device. 

A  comparison  between  Fig.  27  (constant  mobility)  and  Fig.  28 

(field  dependent  mobility)  shows  that  velocity  saturation  produces 

an  Increase  of  saturation  resistance  and,  in  addition,  a  situation 

where  current  saturation  exists  at  a  substantially  lower  drain 

voltage  than  can  be  predicted  by  traditional  theoretical  concepts; 

this  observation  is  qualitatively  consistent  with  laboratory  measure- 
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ments  upon  very  short-channel  IGFET  structures. 

2.4  The  Depletion  Layer  In  an  IGFET 

An  IGFET  contains  three  regions  of  mobile  carrier  depletion: 

1.  the  source  junction  space-charge  layer, 

2.  the  drain  junction  space-charge  layer,  and 

3.  a  portion  of  the  gate  electrode  space-charge  layer. 

These  regions  of  carrier  depletion  have  an  important  function.  Specif¬ 
ically,  they  provide  a  layer  of  semiconductor  material  that  electrically 
insulates  the  conducting  channel  from  a  highly  conductive  substrate.  In 
addition,  the  electrostatic  charge  produced  by  mobile  carrier  depletion 
significantly  influences  the  potential  distribution  in  an  IGFET.  Al¬ 
though  elementary  theory  treats  on  a  one-dimensional  basis  all  depletion 
layer  calculations,  this  is  clearly  an  approximation.  The  geometrical 
configuration  of  an  IGFET  produces  regions  of  depletion  layer  inter¬ 
action  that  must  be  calculated  using  a  two-dimensional  mathematical 
model;  some  initial  investigations  into  this  topic  are  presented  here. 

Figure  29  illustrates  the  calculated  depletion  layer  boundary  in 
a  5.0  um  IGFET  with  an  assumed  gate  bias  of  2.0  volts  and  an  assumed 
drain  bias  of  5.0  volts.  Throughout  this  calculation  the  depletion 
layer  edge  is  defined  as  that  location  where  the  mobile  carrier  den¬ 
sity  is  one-half  the  substrate  impurity  atom  density.  The  overall 
shape  of  this  depletion  layer  edge  is  a  consequence  of  space-charge 
interaction;  two  depletion  layers  (for  example  the  gate  and  drain) 
are  located  within  the  same  region  of  the  substrate  material. 

A  one -dimensional  mathematical  model  adequately  characterizes 
the  depletion  layers  in  regions  of  this  structute  containing  no 
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space-charge  interaction..  In  contrast,  space-charge  layer  interaction 
in  an  IGFET  is  a  consequence  of  two  depletion  layers,  at  right  angles 
to  each  other,  sharing  the  ionized  impurity  atoms  within  a  given  region 
of  the  substrate.  One  consequence  of  this  situation  is  a  substantial 
reduction  in  the  number  of  ionized  impurity  atoms  available  to  each 
depletion  layer;  this  situation  is  equivalent  to  a. reduction  of  sub¬ 
strate  doping.  An  important  consequence  of  space-charge  layer  inter¬ 
action  is  an  increase  of  depletion  layer  width,  and  a  modification  of 
the  IGFET  potential  distribution. 

Before  discussing  the  two-dimensional  aspects  of  these  depletion 
lexers  we  shall  consider  the  one-dimensional  space-charge  layer  width 
due  to  the  gate  electrode,  Fig.  30.  For  this  calculation  the  source- 
drain  voltage  is  zero  and  the  depletion  layer  width  is  established  for 
a  region  far  removed  from  either  the  source  or  drain  junction.  As 
before,  the  space-charge  layer  edge  is  defined  as  that  location  where 
the  majority  carrier  density  is  one-half  the  substrate  impurity  atom 
density. 

From  Fig.  30,  there  is  no  region  of  this  curve  where  the  depletion 
layer  width  is  ideally  independent  of  the  applied  gate  voltage.  This 
situation  could  contribute  to  well  known  inconsistencies  in  the 
definition  of  threshold  voltage.  Elementary  theory  of  IGFET  operation 
is  based  upon  an  assumption  that  the  depletion  layer  width  is  constant 
if  vthiV  when  the  gate  voltage  exceeds  threshold  it  is  assumed 
all  additional  space-charge  arises  from  mobile  carriers  within  the 
source-drain  channel.  Furthermore,  it  is  assumed  that  this  condition 
of  constant  space-charge  width  is  attained  when  the  channel  carrier 
density  equals  the  substrate  Impurity  atom  density.  For  the  5.0  urn 
model  used  in  this  investigation,  threshold  should  occur  at  a  gate  bias 
of  about  0.889  V;  from  Fig.  30  the  depletion  layer  width  at  0.889  volts 
is  only  90%  of  its  maximum  value,  fhis  situation  could  represent  one 
source  of  difficulty  in  defining  the  threshold  voltage. 

In  a  two-dimensional  model,  the  shape  of  interacting  depletion  lexers 
can  be  determined  from  the  direction  and  magnitude  of  the  space-charge 
<*yer  electric  field.  For  example,  regions  of  space-charge  layer 
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Fig.  30.  Calculated  gate  electrode  space-charge  layer  width  In 
a  10.0  vm  MOSFET  when  V$D  =  0. 
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interaction  between  the  gate  and  drain  exhibit  two  electric  field  ■ 
components:  one  parallel  and  one  perpendicular  to  the  oxide  and 
semiconductor  interface.  The  electric  field  .component  parallel  to 
this  interface  is  attributable  to  the  drain  junction^  and  the  dep¬ 
letion  layer  edge  is  assumed  to  be  located  where  this  electric  field 
component  has  decreased  to  1000  V/ cm.  Figure  31  illustrates  this 
calculated  depletion  layer  edge  when  VSG  =  2.0  volts  and  VSD  =  0.  > 

Figures  32-35  show  the  calculated  drain  junction  space-charge' 
layer  edge  with  an  increase  of  source-drain  biasing  voltage.  This 
sequence  of  calculations  exhibit  an  expansion  of  the  drain  junction 
depletion  layer  into  the  gate  electrode  depletion  layer.  As  previously 
stated,  this  expansion  is  attributed  to  the  lack  of  ionized  impurity 
atoms  available  to  the  drain  junction;  a  large  percentage  of  these 
ions  are  electrostatically  tied  to  the  gate  electrode.  For  this  reasor), 
this  gate  depletion  layer  appears  as  a  region  of  small  impurity  atom 
density  through  which  the  drain  depletion  laydr  can  reach  through  to 
the  source  junction. 

From  the  suggested  mechanism  producing  reach-through,  the' degree  ■ 
of  drain  depletion  layer  expansion  into  the  gate  depletion  layer  should 
be  dependent  upon  the  source-gate  voltage;  this  is  confirmed  by  the 
calculations  shown  in  Figs.  36-38.  In  this  sequence  of  calculations 
(Figs.  36-38)  the  source-drain  bias  is  maintained  at  5.0  volts  and  the 
source-gate  voltage  has  a  magnitude  of  1.0,  2.0,  and  4.0  volts,  res-1 
pectively.  At  a  gate  bias  of  1.0  volt  (Fig.  36)  there  is  little  penet¬ 
ration  of  the  drain  space-charge  layer  into  the  gate  space-charge  layer. 
In  contrast,  at  a  source-gate  bias  of  4.0  volts  (Fig.  3.8)  the  drain 
depletion  layer  penetrates  to  the  source  junction  and  produces  a 
contraction  of  the  source  junction  space-charge  layer. 

At  this  time  no  compandor,  has  been  made  between  the  reach-through 
mechanism  suggested  by  Figs.  32-35  and  traditional  source-drain 
punch- through.  In  traditional  punch-through  the  drain  space-charge 
layer  penetrates  to  the  source  and  produces  a  contraction  of  the  source 
junction  space-charge  layer;  this  contraction  is  viewed  as  equivalent 
to  the  application  of  a  forward  bias  to  the  source-junction.  Figure  37 
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shows  a  comparison  between  the  drain  junction  space-charge  1«yer  edge 
when  VSD  =  0  and  also  when  VSD  =  5.0  volts.  Clearly,  from  Fig.  37, 
reach-through  produces  a  contraction  of  the  source  junction  space-charge 
layer,  although  the  influence  of  reach- through  upon  the  source-drain 
current  is  a  topic  for  further  investigation. 


2.5  The  Threshold  Voltage  in  an  I6FET 

18 

From  elementary  one-dimensional  theory,  the  source-drain  electric 
current  in  the  low  current  region  of  an  IGFET  volt-ampere  characteristics 
is  given  by 


*SD  =  ^^SG^th^SD*  1  VSD^ 


where 


n  = 


Weoxu 

Lctox 


(7) 

(8) 


In  the  saturation  region  of  these  volt-ampere  characteristics  this 
source-drain  electric  current  is  given  by 


where 
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and 
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Figure  39  presents  a  comparison  between  this  elementary  IGFET 
theory  (eqs.  7  and  9)  and  a  rigorous  calculation  of  these  volt-ampere 
curves.  At  each  gate  voltage  the  threshold  voltage  (V^  in^gs^JL-and- 
9)  was  selected  to  prodirce^greemeTTtn^tfi-ffi^^  drain 

current  at  a  drain  biasing  voltage  producing  deep  saturation.  From 
Fig.  39,  substantial  agreement  is  obtained  between  these  two  theoretical 
techniques,  although  to  attain  this  agreement  a  different  threshold 
voltage  was  required  at  each  value  of  applied  gate  voltage. 


From  Fig.  39  the  inferred  threshold  voUage  changes  from  1.09  volts 
to  1.21  volts,  depending  upon  V<.g.  The  reasons  for  this  situation  are 
presently  unknown;  it  could  result  from  approximations  used  in  the 
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development  of  this  elementary  theory  or,  instead,  from  the  fact  that 
this  elementary  theory  remains  incomplete.  It  appears  that  present 
methods  for  determining  Vth  are  inconsistent  and  the  topic  requires 
further  investigation. 

2.6  The  Potential  Distribution  in  an  IGFET 

In  the  analysis  of  IGFET  operation,  the  potential  distribution 
along  a  source-drain  channel  has  been  subject  to  extensive  investigation. 
Most  workers  approach  this  problem  on  a  one-dimensional  basis,  using  a 
gradual  channel  type  of  approximation  that  was  initially  developed  by 
Shockley  for  the  analysis  of  JFET  structures.  This  one-dimensional 
analysis  has  been  assumed  adequate  (without  extensive  verification)  for 
structures  biased  well  below  current  saturation,  and  for  channel  poten¬ 
tial  distribution  calculations  between  the  source  junction  and  the 
pinch-off  point  (when  biased  into  saturation).  The  two-dimensional 
aspects  of  this  problem  are  recognized  to  be  important  in  calculations 
of  the  potential  distribution  between  the  pinch-off  point  and  the  drain, 
although  there  is  no  simple  method  whereby  this  potential  can  be 
established. 

The  difficulties  associated  with  solving  this  two-dimensional 
problem  are  evident  in  the  publication  of  Chiu  and  Sah.4  Specifically, 
Laplace's  equation  was  solved  for  the  region  of  insulating  material 
that  is  approximately  bounded  by  the  channel  pinch-off  point  and  the 
drain.  Next,  this  solution  equation  was  matched  to  a  solution  of 
Poisson's  equation  for  the  potential  distribution  adjacent  the  insulator 
semiconductor  boundary.  In  this  manner,  a  two-dimensional  analytical 
expression  was  developed  for  the  potential  distribution  in  the  drain 
region  of  an  IGFET. 

It  is  emphasized,  this  mathematical  method  produces  a  very  difficult 
problem  of  matching  (in  magnitude  and  derivative)  two  different  (and 
complicated)  solution  equations  at  the  insulator-silicon  interface.  A 
consequence  of  this  situation  is  the  necessity  of  introducing  approx¬ 
imations  that  simplify  this  mathematical  problem  without  decreasing  the 
generality  of  the  solution.  In  the  work  of  Chiu  and  Sah,  an  approx¬ 
imation.  inadvertently  introduced  that  removed  from  their  solution  \ 
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equation  the  influence  of  ionized  impurity  atoms  upon  the  calculated 
potential  along  the  insulator-semiconductor  boundary.  The  consequences 
of  this  approximation  method  are  yet  to  be  determined. 

Figures  40  and  41  illustrate  computer  calculated  two-dimensional 
potential  distributions  within  a  5.0  ym  I6FET.  Figure  40  shows  this 
potential  distribution  when  the  structure  is  biased  below  electric 
current  saturation;  in  Fig.  41  the  structure  is  assumed  biased  well 
into  current  saturation.  The  two-dimensional  nature  of  these  potential 
distributions  is  a  consequence  of  electrostatic  interactions  of  both 
the  gate  electrode  and  the  drain  in  association  with  fixed  lattice 
charges  within  a  given  region  of  the  substrate.  Regions  of  the  drain 
junction  far  removed  from  the  gate  electrode  exhibit  a  potential 
distribution  that  is  accurately  characterized  by  one-dimensional  abrupt 
p-n  junction  theory.  In  contrast,  near  the  insulator  and  semiconductor 
interface  fixed  lattice  charges  are  shared  between  the  gate  and  drain; 
therefore,  the  contours  of  constant  potential  are  neither  parallel  to 
the  gate  electrode  nor  the  drain  junction. 

Substantial  differences  can  be  observed  between  the  potential  dis¬ 
tribution  shown  in  Figs.  40  and  41.  In  particular.  Fig.  41  shows  a 
point  along  the  semiconductor  to  oxide  interface  at  which  the  semicon¬ 
ductor  potential  just  equals  the  gate  potential;  on  the  source  side  of 
this  location  the  electric  field  normal  to  the  interface  will  hold 
mobile  carriers  against  tne  semiconductor  surface  (thereby  forming  the 
channel)  and  on  the  drain  side  this  field  forces  carriers  away  from  the 
surface.  It  is  emphasized,  within  both  the  conductive  channel  and  the 
depletion  layer  adjacent  this  channel  a  substantial  electric  field  can 
be  shown  to  force  these  mobile  carriers  toward  the  drain  junction. 

In  the  source  portion  of  the  source-drain  conductive  channel 
(Fig.  41)  the  electric  field  at  ilu:  si 1 icon-oxide  interface  forces 
mobile  carriers  against  the  semiconductor  surface;  this  field  is 
attributed  to  an  electrostatic  interaction  between  the  gate  junction 
and  electrostatic  charges  within  both  the  channel  (mobile  carriers) 
and  the  depletion  layer  (ionized  impurity  atoms).  In  the  drain  portion 
of  the  source-drain  channel  this  electric  field  is  in  the  opposite 
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direction;  it  forces  mobile  carriers  away  from  the  semiconductor 
surface.  Clearly,  this  reversal  of  electric  field  results  from  an 
electrostatic  polarity  reversal  between  the  gate  electrode ,and  the  i  ' 
semiconductor  material.  As  a  consequence,  in  this  drain  region, the 
electrostatic  charges  seen  by  the  gate  electrode  must  be  of  opposite 
polarity  to  those  seen  in  the  source  region,  of  this  channel.  Figure 
42  shows  how  this  is  accomplished. 

Figure  42  illustrates  the  calculated  contours  of  constant  relative 
electrostatic  displacement  between  the  gate  electrode  and  drain.  The 
electric  field  attributable  to  this  mechanism,  in  conjunction  with  the 
electric  field  produced  by  ionized  impurities  in  the  semiconductor 
material,  produce  the  potential  distribution  shown  in  Fig.  41  (at  the 
drain  end  of  the  channel).  From  Fig.  42  some  charges  within  the  drain 
region  are  electrostatically  tied  to  the  gate  electrdde,  ,and  thereby.we 
obtain  the  required  reversal  of  electric  field. 

From  Figs.  40  and  41,  it  becomes  evident  that  several  weaknesses 
should  be  observed  in  a  one-dimensional  approximation  for  the  potential 
distribution  along  a  source-drain  channel.  For  example,  at  large  values 
of  drain  voltage  this  approximate  theory  becomes  questionable:  at 
voltages  sufficient  to  produce  substantial  drain  junction  space-charge 
layer  expansion  into  the  gate  depletion  layer.  Similarly,  between  the 
point  of  channel  termination  and  the  drain  junction  one-dimensional 
theory  takes  into  consideration  only  the  potential  distribution  arising 
from  the  drain  junction  space-charge  layer.  These  weaknesses  of  the 
approximate  theory  become  evident  in  Figs.  43-48. 

! 

Figures  43-45  present  a  comparison  between  one-dimensional  and 
two-dimensional  calculations  of  source-drain  channel  potential 
distribution  for  the  10.0,  5.0  and  1.0  urn  structure,  respectively. 

From  these  illustrations,  subshaiu  •  1 'y  better  agreement  is  obtained  at 
small  values  of  source-drain  voltage.  It  is  suspected  that  this  agree¬ 
ment  could  be  significantly  improved  if  the  source  junction  space-charge 
layer  potential  distribution  was  included  in  the  one-dimensional  approx¬ 
imation.  At  large  values  of  drain  voltage  the  one-dimensional  theory  is 
shown  to  be  incomplete  in  its  characterization  of  the  channel  potential 
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Calculated  contours  of  constant  relative  electric 
displacement  between  the  gate  electrode  and  the  drain. 
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distribution.  At  best,  this  one-dimensional  theory  provides  only  a 
first-order  type  of  approximation,  and  should  be  used  with  caution 
in  the  evaluation  of  IGFET  structures. 

Similarly,  Figs.  46-48  present  a  comparison  between  these 
one-dimensional  and  two-dimensional  calculations  when  the  source- 
drain  voltage  is  held  constant  (VSD  =  3.0  volts)  and  the  source-gate 
voltage  is  changed  from  1.0  to  4.0  volts.  At  small  values  of  gate 
voltage,  the  appreciable  error  in  this  approximate  theory  appears  to 
be  attributable  to  a  poor  characterization  of  the  channel  potential 
distribution  near  the  drain  end  of  the  structure.  As  previously  shown, 
an  increase  of  gate  voltage  produces  a  significant  increase  of  drain 
junction  penetration  along  the  source-drain  channel;  a  situation  that 
has  not  been  approximated  in  this  one-dimensional  theory. 
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CHAPTER  III 

Bipolar  Transistor  Investigations 
3.0  Introduction 

Important  questions  remain  concerning  the  physical  mechanisms 
involved  in  bipolar  transistor  operation.  In  particular,  there  is 
extensive  technical  literature  on  the  theory  of  current  gain  in 
bipolar  transistor  operation  (both  high  frequency  and  low  frequency), 
although  a  detailed  comparison  with  experiment  shows  this  theory  has 
many  deficiencies.  Most  comparisons  between  theory  and  experiment 
show  the  measured  ft  of  a  transistor  is  two-thirds  to  one-half  its 
calculated  value.  This  difference  between  experiment  and  theory 
appears  to  arise  from  an  inherent  inadequacy  of  our  mathematical 
models:  the  theory  of  bipolar  transistors  is  based  upon  one-dimen¬ 
sional  and  pseudo  two-oimensional  models  that  do  not  take  into 
consideration  detailed  two-dimensional  mechanisms  encountered  in 
this  semiconductor  device. 

Numerous  workers  have  conducted  a  pseudo  two-dimensional  analysis 
of  bipolar  transistor  operation  by  coupling  together  several  one-dimen¬ 
sional  solutions  through  either  a  base  region  resistor^  or  through  more 

2 

complicated  base  region  mechanisms.  This  technique  is  an  obvious 
attempt  to  overcome  computational  difficulties  inherent  in  any  truly 
two-dimensional  numerical  solution  of  this  problem.  Between  six  and 
twelve  one-dimensional  arrays  are  used  to  describe  the  intrinsic  region 
of  a  bipolar  transistor,  and  these  arrays  are  coupled  together  at  three 
locations:  in  the  emitter  region,  in  the  collector  region,  and  through 
a  base  region  element.  This  technique  has  the  advantage  of  computational 
simplicity,  yet  there  is  little  information  available  to  establish  the 
errors  introduced  by  such  an  approximation  method. 

It  is  for  these  reasons  we  have  undertaken  the  task  of  developing 
a  rigorous  two-dimensional  solution  of  the  bipolar  transistor  problem. 

The  complications  associated  with  running  this  program  are  recognized 
to  render  it  inapplicable  for  the  solution  of  engineering  problems  on 
a  routine  basis.  The  intent  here  is -to  use  this  program  in  an  inves¬ 
tigation  of  the  mechanisms  involved  in  transistor  operation  and,  thereby. 
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obtain  a  verification  of  the  more  approximate  computational  techniques 
presently  in  use. 

Our  computer  program  utilizes  a  matrix  of  2800  points  that  are 
spatially  distributed  to  reduce  errors  associated  with  approximating 
a  transistor  in  this  fashion.  An  important  difficulty  is  the  com¬ 
putation  time  required  to  attain  a  single  solution;  this  computation 
time  varies  significantly  with  the  assumed  operating  level  of  the 
transistor.  At  low  levels  of  emitter  current  only  a  small  amount  of 
base  region  conductivity  modulation  is  realized  and,  as  a  consequence, 
a  solution  is  attained  in  a  few  minutes  (360/Mod  91).  In  contrast, 
at  large  values  of  emitter  current  the  structure  exhibits  substantial 
conductivity  modulation;  this  situation  can  increase  the  computation 
time  to  in  excess  of  4.0  hours. 

Although  the  computation  time  is  excessive,  using  this  bipolar 
transistor  program,  a  long  computation  time  is  not  inconsistent  with 
previous  experience.  All  programs  developed  during  the  past  few  years 
undergo  several  phases  of  development;  the  initial  prototype  always 
requires  an  excessive  amount  of  computation  time  and,  after  further 
program  development,  the  required  computation  time  decreases  signif¬ 
icantly.  For  this  reason,  before  we  can  use  this  computer  program 
for  an  extensive  investigation  of  bipolar  transistor  operation,  further 
work  is  necessary  to  reduce  the  computation  time  to  an  economically 
practical  level. 

The  present  discussion  is  intended  to  outline  qualitative  obser¬ 
vations  based  upon  initial  test  calculations  using  this  computer 
program  One  sequence  of  calculations  was  conducted  for  a  typical 

bipolar  transistor  structure,  throughout  a  range  of  emitter  current 

.  .  4  2 

densities  up  to  about  3.8  x  10  amps /cm  .  The  principal  purpose  for 

this  effort  was  to  establish  the  operating  characteristics  of  the 

program  and,  as  a  secondary  issue,  to  qualitatively  study  some  concepts 

of  transistor  operation  that  are  in  question.  In  particular,  an  attempt 

was  made  to  qualitatively  establish  the  degree  of  emitter  current  ■ 

crowding  obtained  in  a  bipolar  transistor  at  large  emitter  current 

densities.  In  conjunction  with  these  calculations,  a  qualitative  study 

was  also  made  to  establish  the  mechanisms  of  current  gain  fail-off  at 

large  values  of  emitter  current  density. 
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3.1  Analysis 

It  has  been  shown  that  the  hole  and  electron  distributions  in 

q 

semiconductor  material  are  described  by  the  equations 


a)  div  grad  if>  =  ^3-  (N  -  n  +  p) 

K£q 

b)  Jp  =  -  qDp  grad  p  -  qPpp  grad  * 

c)  ?n  *  qDn  grad  n  -  qpnn  grad  (1) 


d) 

e) 

f) 


div  Jp  =  qRp 

div  Jn  =  qR„ 

1  ■  1  +  J„ 
T  p  n 


when  it  is  assumed  that  no  trapping  mechanisms  exist  within  the 
structures  under  consideration.  Finite  difference  methods  are  used 
co  numerically  solve  this  system  of  equations;  details  of  these 
computational  methods  are  outlined  in  Chapter  IV. 


As  in  all  mathematical  investigations  of  bipolar  transistor 
operation,  little  information  is  available  concerning  the  minority 
carrier  lifetimes  within  the  base  and  emitter  region.  Therefore, 
for  this  initial  sequence  of  calculations  we  have  assumed  the  hole 
(emitter)  and  electron  (base)  minority  lifetimes  are  identical, 
and  the  magnitude  of  these  lifetimes  was  adjusted  to  yield  a  cal¬ 
culated  current  gain  of  about  200;  this  was  accomplished  with  an 
assumed  lifetime  of  10  nsec. 


Carrier  mobilities  in  these  calculations  are  dependent  upon 
local  impurity  atom  densities,  local  electric  field,  and  direction.^ 
The  characterization  of  this  physical  parameter  at  large  values  of 

v 

electric  field  is  consistent  with  published  information  on  this 
5 

topic.  Specific  information  on  the  methods  used  to  describe  this 
minority  carrier  mobility  is  outlined  in  Chapter  IV. 

Throughout  this  mathematical  investigation,  particular  emphasis 
is  placed  upon  the  selection  of  boundary ^conditions  that  do  not  intro 
duce  errors  in  the  calculated  results.  In  particular,  the  boundary 
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conditions  used  in  this  analysis  approximate  physical  and  electrical 
characteristics  at  the  outer  periphery  of  this  semiconductor  structure, 
rather  than  along  internal  boundaries  established  by  its  physical  or 
electrical  properties.  The  exposed  semiconductor  surface  is  assumed 
to  be  an  ideal  electrical  insulator;  i.e.,  no  electric  current  is 
permitted  normal  to  these  bounding  surfaces.  The  ohmic  contacts  are 
approximated  by  equi potential  surfaces  that  are  charge  neutral;  where 
appropriate,  these  ohmic  contacts  have  been  located  sufficiently  far 
from  the  active  regions  of  this  device  to  have  no  influence  upon  its 
electrical  properties. 

Figure  49  shows  the  analytical  model  used  in  these  two-dimensional 
transistor  calculations.  In  this  model  we  have  eliminated  a  large  por¬ 
tion  of  the  extrinsic  base  region  and,  thereby,  concentrated  the  available 
matrix  points  in  the  intrinsic  portion  of  the  structure.  It  is  emphasized, 
the  base  ohmic  contact  location  in  Fig.  49  is  a  sufficient  distance  from 
the  emitter  junction  so  that  this  contact  has  no  influence  upon  the 
fundamental  mechanisms  of  transistor  operation. 

The  base  region  in  this  structure  is  assumed  to  result  from  a 
two-step  diffusion6  into  n-type  semiconductor  material  of  homogeneous 
background  doping  (10  atoms/cm  ).  For  this  model  the  boron  profile 

is  a  one-dimensional  Gaussian  distribution  with  a  surface  concentration 

19  3 

of  10  atoms/cm  .  The  emitter  distribution  has  been  calculated  using 
a  5  pm  diffusion  mask  opening;  the  profile  calculation  is  based  upon 
a  previous  publication  on  this  topic. ^  It  is  assumed  this  emitter 
junction  arises  from  a  phosphorus  surface  concentration  of  5  x  10  atoms/cm 
that  is  diffused  for  a  sufficient  period  to  produce  an  emitter  junction 
depth  of  1.5  urn.  All  calculations  are  based  upon  an  assumed  emitter 
junction  length  of  10.0  pm. 


3.2  Current  Gain  in  a  Bipolar  Transistor 


At  relatively  small  values  of  emitter  current,  it  is  experimentally 
observed  that  the  current  gain  of  a  transistor  increases  with  an  increase 
of  emitter  current.  In  contrast,  at  large  values  of  current  this  gain  is 
observed  to  decrease  with  an  increase  of  emitter  current.  These  experi¬ 
mental  observations  are  consistent  with  our  calculations  of  current  gain 
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Fig.  49.  Analytical  model  for  a  bipolar  transistor 


(Fig.  50),  although  it  will  be  shown  that  traditional  explanations 
for  this  situation  are  Inconsistent  with  conclusions  drawn  from  our 
two-dimensional  analysis. 

In  the  low  current  region  of  Fig.  50,  an  Increase  of  current  gain 
(with  an  Increase  of  emitter  current)  Is  traditionally  attributed  to 
ambl polar  diffusion  of  minority  carriers  across  a  transistor  base 

O 

region.  This  mechanism  was  first  proposed  by  W.  M.  Webster  (In  1954) 
and  has  been  generally  accepted  throughout  the  semiconductor  industry. 
Briefly,  It  Is  presumed  that  conductivity  modulation  within  the  tran¬ 
sistor  base  region  produces  an  electric  field  that  enhances  minority 
carrier  transport  from  emitter  to  collector.  At  large  values  of 
conductivity  modulation,  It  can  be  shown  this  electric  field  enhance¬ 
ment  Is  equivalent  to  doubling  the  base  region  minority  carrier  diffusion 
coefficient.  An  Increase  of  base  region  transport  efficiency  (due  to 
ambi polar  diffusion  mechanisms)  Is  presumed  to  produce  the  observed 
increase  of  transistor  gain,  with  an  Increase  of  emitter  current. 

Although  ambi polar  diffusion  mechanisms  do  take  place  In  the 
transistor  base  region,  rigorous  calculations  show  that  this  mechanism 
does  not  necessarily  produce  the  gain  vs.  emitter  current  characteris¬ 
tics  observed  in  diffused  transistor  operation.  Figure  51  illustrates 
the  one-dimensional  impurity  atom  distribution  in  a  typical  diffused 
bipolar  transistor,  included  In  this  illustration  are  arrows  (marked  J) 
indicating  the  component  of  minority  carrier  flux  (electrons)  attrlbutabl 
to  drift  mechanisms  within  the  transistor  base  region,  at  small  levels 
of  base  region  conductivity  modulation.  These  arrows  show  the  base 
region  electric  field  In  a  diffused  transistor  does  not  necessarily 
enhance  minority  carrier  transport  between  the  emitter  and  collector, 

q 

as  presumed  in  the  theory  of  drift  transistors.  Near  the  emitter 
junction,  compensation  of  base  region  impurity  atoms  by  emitter  region 
impurity  atoms  yields  a  "built-in"  electric  field  that  opposes  minority 
carrier  transport  from  emitter  to  collector;  this  opposition  produces 
a  reduction  of  base  region  transport  efficiency. 

A  relatively  small  Impurity  atom  density  exists  throughout  this 
region  of  reverse  "built-in"  electric  field;  therefore,  this  region 
is  highly  susceptible  to  conductivity  modulation.  Calculations  indicate 
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Fig.  51.  Illustrative  one-dimer 
for  a  diffused  bioola> 
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that  at  small  values  of  emitter  current  conductivity  modulation 
essentially  eliminates  this  reverse  ’’built-in"  field  and,  therefore, 
increases  the  base  region  transport  efficiency.  Thereby,  an  increase 
of  transistor  gain,  with  an  increase  of  emitter  current,  is  obtained 
at  relatively  small  values  of  emitter  current. 

Traditional  theory  predicts  that  at  large  emitter  current  den¬ 
sities  bipolar  transistors  will  exhibit  a  decrease,  of  current  gain, 
with  an  increase  of  emitter  current,  although  this  theory  is  based 
upon  a  one-dimensional  analysis  of  the  problem.  Our  calculations  of 
current  gain  (Fig.  50)  are  consistent  with  this  prediction,  yet  the 
mechanisms  producing  a  decrease  of  gain  are  essentially  two-dimensional. 
From  these  calculations,  it  is  suggested  that  a  substantial  modification 
is  necessary  In  the  traditional  theory  of  bipolar  transistor  operation 
before  this  theory  is  entirely  consistent  with  the  physical  mechanisms 
Involved  in  device  operation.  , 

Figure  52  shows  the  calculated  base  region  transport  efficiency 
(6g)  and  emitter  region  injection  efficiency  ( e^)  for  this  analytical 
model.  These  parameters  ( eB ,  6^)  represent  the  current  gain  that  would 
be  exhibited  by  this  transistor  if  the  individual  parameter  under  con¬ 
sideration  represented  the  Only  source  of  gain  lin’tation,  For  example, 
3g  represents  the  current  gain  if  the  transistor  contained  ar.  ideal 
emitter  junction  (y=1).  Similarly,  6^  represents  the  transistor  current 
gain  assuming  no  minority  carrier  recombination  mechanisms  within  the 
base  region.  1 

From  Fig.  52,  botn  emitter  injection  efficiency  (e^)  and  base 
region  transport  efficiency  (6g)  undergo  a  substantial  decrease  at 
large  values  of  emitter  current.  Traditional  theory  would  attribute 
this  loss  of  base  transport  efficiency  ( 6g)  to  the  so-called 
"Kirk-effect":  collector  junction  displacement  (and,  hence,  base 
region  widening)  due  to  a  large  minority  carriep  density  within  the 
collector  junction  space-charge  layer. In  addition,  this  traditional 
theory  would  attribute  the  loss  of  injection  efficiency  (e^)  to  base 
region  conductivity  modulation.  Although  any  two-dimensional  structure 
must,  indeed,  exhibit  the  mechanisms  suggested  by  this  traditional 
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Fig.  52.  Calculated  emitter  transport  and  base  region  transport 
efficiencies  for  a  bipolar  transistor  (Fig.  49). 
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theory,  there  Is  no  reason  to  assume  that  these  mechanisms  produce 
the  observed  current  gain  characteristics  of  a  diffused  bipolar 
transistor. 

Calculations  indicate  that  base  region  widening,  due  to  minority 
carriers  within  the  collector  space-charge  layer  (Kirk -effect) ,  is  of 
little  consequence  In  this  particular  semiconductor  device.  Instead, 
the  calculated  decrease  of  base  region  transport  efficiency  (eB) 
appears  to  result  from  side-wall  Injection  of  minority  carriers  by 
the  emitter  junction.  This  situation  is  Illustrated  in  Fig.  53. 

In  Fig.  53,  a  dashed  line  Is  used  to  designate  that  region  of  the 
emitter  junction  where  the  current  density  (due  to  electron  injection) 
is  greater  than  one-half  its  maximum  value.  At  small  values  of  total 
emitter  current  (0.19  ma)  only  a  small  amount  of  electron  injection 
takes  place  along  the  emitter  side-wall.  In  contrast,  at  large  values 
of  total  emitter  current  (24  ma)  substantial  side-wall  injection  Is 
encountered  In  this  semiconductor  device. 

A  consequence  of  this  side-wall  Injection  is  an  Increase  of 

effective  base  width,  with  an  increase  of  emitter  current  (see  Fig.  54). 

Carriers  entering  the  base  region  from  the  emitter  junction  side-wall 

must  traverse  a  longer  path  than  carriers  entering  from  other  regions 

of  the  emitter;  these  side-wall  injected  carriers  exhibit  substantial 

base  region  recombination  and,  thereby,  produce  the  calculated  fall -off 

of  base  transport  efficiency  (eg).  It  is  emphasized,  at  high  emitter 

currents  the  region  of  base  region  recombination  differs  significantly 

from  that  suggested  by  traditional  theory  of  bipolar  transistors. 

Traditional  theory  suggests  that  most  base  region  recombination  takes 

place  under  the  emitter,  rather  than  along  its  periphery,  and  this 

concept  has  been  adopted  in  most  pseudo  twc-dimensional  investigations 

1  2 

of  device  operation. 

As  previously  stated,  it  is  generally  assumed  that  the  decrease 
of  emitter  injection  efficiency  with  emitter  current  (e^  in  Fig.  52) 
is  a  result  of  base  region  conductivity  modulation.  All  one-dimensional 
(and  pseudo  two-dimensional)  calculations  of  device  performance  are 
based  upon  an  assumption  that  this  is  the  dominant  mechanism  producing 
&  fall -off.  Although  base  region  conductivity  modulation  may,  indeed, 


87 


Depth  ifim) 


Senmnductor  Surface 


Distance  from  Device  Center  ifim) 


Pig,  53.  Illustration  of  the  calculated  active  emitter  junction 
area  with  increasing  emitter  current 


m 


Low-Current  Injection 


Collector 


Collector 

Fig.  54.  Illustration  of  the  calculated  minority  carrier  flux 
distribution  in  a  bipolar  transistor. 
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produce  a  loss  of  emitter  injection  eff icier cy,  Fig.  55  shows  another 
mechanism  that  has  a  more  significant  influence  upon  the  electrical 
characteristics  of  this  particular  bipolar  transistor. 

Figure  55  shows  a  sketch  of  this  two-dimensional  emitter  junction 
upon  which  is  drawn  dashed  lines  to  indicate  the  region  of  maximum  back 
injection  (holes)  into  the  emitter.  At  each  value  of  total  emitter 
current,  a  given  dashed  line  bounds  that  portion  of  the  junction  exhibit¬ 
ing  a  back  injection  (holes  Injected  into  the  emitter  region)  greater 
than  one-half  its  maximum  value.  At  a  small  value  of  emitter  curreit, 
little  back  injection  takes  place  at  the  emitter  side-wall.  In  contrast, 
at  large  values  of  emitter  current  (14  ma)  almost  all  the  loss  of  Injec¬ 
tion  efficiency  is  attributable  to  side-wall  Injection  of  holes  into  the 
emitter  region. 


A  consequence  of  this  situation  is  that  the  loss  of  emitter  injec¬ 
tion  efficiency  is  attributable  to  side-wall  injection  into  the  emitter 
region,  rather  than  from  the  intrinsic  region  of  this  structure  (Fig.  56). 
It  is  emphasized,  at  large  values  of  emitter  current  the  distribution 
of  base  region  current,  due  to  back  Injection,  differs  significantly  from 
the  distribution  suggested  by  elementary  theory.  At  large  values  of 
emitter  current  the  base  region  flux,  due  to  this  mechanism,  is  concen¬ 
trated  near  the  emitter  periphery,  rather  than  between  the  emitter  and 
collector  junctions. 


The  foregoing  represent  some  Initial  conclusions  drawn  from  this 
two-dimensional  analysis.  These  calculations  show  the  side-wall  of  an 
emitter  junction  has  an  important  influence  upon  the  electrical 
properties  of  this  semiconductor  device.  Further  study  of  the  computed 
results  is  required  before  we  can  quantitatively  establish  the 
mechanisms  producing  side-wall  injection  and,  thereafter,  generate  «n 
elementary  (or  engineering)  theory  that  is  consistent  with  the  physics 
of  bipolar  transistor  operation. 


Depth  (fim) 


Distance  from  Device  Center  (/xm) 


Fig.  55.  Illustration  of  the  calculated  region  of  back  Injection 
by  the  emitter  junction  with  increasing  emitter  current 
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Fig.  56.  Illustration  of  the 
diffused  emitter  in 
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calculated  back  injection  of  a 
a  bipolar  transistor. 


3. 3  Emitter  Current  Crowdin 


Traditional  theory  of  transistor  operation  is  based  upon  an  as¬ 
sumption  that  all  base  region  majority  carrier  current  is  attributable 
to  mechanisms  encountered  between  the  emitter  and  collector  junctions 
(base  region  recombination  and  back  injection  by  emitter  junction),  as 
illustrated  in  Mg.  57.  An  important  consequence  of  this  assumption 
is  a  longitudinal  variation  of  emitter  biasing  voltage  that  produces  a 
substantial  increase-  of  emitter  current  density  near  the  emitter  edge 
(current  crowding).  It  is  shown  here  that  emitter  current  crowding  is 
of  minor  consequence  in  this  particular  mathematical  model;  present 
pseudo  two-dimensional  analysis  of  this  same  structure  indicates  sub¬ 
stantial  current  crowding. 

Figure  58  shows  the  calculated  emitter  crowding  in  this  semicon¬ 
ductor  device,  as  a  function  of  total  emitter  current.  From  this 
sequence  of  calculations,  at  a  total  emitter  current  of  ,13  ma,  this 
structure  exhibits  a  current  crowding  of  only  0.8:  ratio  of  current 
density  at  center  of  emitter  to  current  density  at  emitter  edge.  In 
contrast,  a  pseudo  two-dimensional  calculation  would  show  a  substan¬ 
tially  greater  current  crowding  ratio  at  this  same  total  emitter 
current.  The  reasons  for  this  disagreement  are  e/'dent  in  view  of 
the  current  gain  mechanisms  outlined  in  section  3.3.  ' 

Figure  59  illustrates  the  path  of  majority  carrier  flux  in  this 
transistor  base  region.  At  small  values  of  emitter  current  this 
majority  carrier  flux  appears  under  the  emitter  junction,  although  its 
magnitude  is  insufficient  to  produce  a  significant  amount  of  emitter 
crowding.  At  large  values  of  emitter  current,  base  region  majority 
carrier  flux  is  directed  predominantly  toward  the  emitter  side-wall 
and  recombine tion  centers  located  outside  the  intrinsic  base  region; 
this  base  region  flux  cannot  influence  the  longitudinal  distribution 
of  emitter  bias  and,  hence,  cannot  produce  current  crowding.  As  a 
consequence,  a  rigorous  two-dimensional  solution  of  the  bipolar 
transistor  problem  shows  substantially  less  emitter  current  crowding 
than  is  suggested  by  the  present  technical  literature. 


■WT! 


I 

Fig.  57.  Illustration  of  the  assumed  majority 

carrier  flux  tn  a  transistor  base  region 
(traditional  theory). 
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Fig.  59. 


Illustration  of  base  region  majority  carrier  flux  in 
a  bipolar  transi s tor. 


3..4  References 

1.  H.  N,  Ghosh,  IEEE  Trans .,  ED-12,  513  (1965) 

2.  B.  V.  Gokhale,  IEEE  Trans.,  ED-17,  594  (1970) 

3.  W.  Van  Rocrbroeck,  Bell  Syst.  Tech.  7.,  29,  560  (1950) 

4.  D.  M.  G.  jhey  and  R.  F.  Thomas,  Proc.  IEEE  (Letters),  55, 

2192  (1967) 

5.  E.  J.  Ryder,  Phys.  Rev.,  90,  766  (1953) 

6.  D.  P.  Kennedy  and  P.  C.  Murley,  Proc.  IEEE,  5'!,  620  (1964) 

7.  D.  P.  Kennedy  and  R.  R.  O'Brien,  IBM  J.  Res.  and  Dev.,  9,  179  (1965) 

8.  W.  M.  Webster,  Proc.  IRE,  42,  914  (1954) 

9.  H.  Kroner,  Naturwissenschaften,  40,  578  (1953) 

10.  C.  T.  Kirk,  IRE  Trans.,  ED-9,  164  (1962) 


\ 


97/98 


CHAPTER  IV 


Computational  Methods 

4.0  Introduction 

In  this  chapter  we  describe  the  two-dimensional  mathematical  methods 
used  throughout  this  investigation.  Our  methods  are  based  on  a  well-known 
system  of  partial  differential  equations  describing  steady-state  carrier 
distributions  in  semiconductors.^  The  geometry  and. boundary  conditions 
adopted  to  approximate  the  structures  under  consideration  are  intended 
to  be  consistent  with  idealized  devices.  Finite  difference  techniques 
are  applied  to  this  system  of  differential  equations;  an  algorithm  is 
presented  for  computing  solutions  for  these  equations  as  applied  to  an 
IGFET  and  bipolar  transistor. 

Since  several  papers  have  appeared  on  the  use  of  difference  schemes 
for  semiconductor  problems,  this  discussion  will  be  limited  to  the 
individual  details  of  the  present  method.  However,  some  of  our  results 

2  3 

do  not  agree  with  previously  published  computations  on  IGFET  structures;  ’ 
we  therefore  present  enough  details  of  the  computations  to  justify  our 
results. 

Our  analytical  model  for  an  IGFET  is  shown  in  Fig.  60.  This 
structure  is  divided  into  two  regions:  first,  the  substrate  ( R-j )  and, 
second,  the  oxide  layer  (R2).  The  electrostatic  potential  function  is 
defined  in  the  region  R=R^UR2;  the  carrier  concentrations,  current 
densities,  etc.,  are  defined  only  in  the  substrate  region  R^. 

A  number  of  approximations  have  been  introduced  into  this  analysis, 
to  simplify  the  co\\  jstions  and  the  interpretation  of  the  computed 
results.  A  n-channe!  a?rce  is  considered,  in  which  hole  current 
and  recombination  are  neale«-ted.^  We  also  assume  abrupt  junctions  at 
the  source  and  drain,  a  uniform  -ondegenerate  substrate  doping,  and  an 
ideal  charge-free  substrate -oxide  interface.  With  these  simplifications, 
the  mathematical  problem  considered  here  is  reduced  to  the  simultaneous 
solution  of  the  Poisson  equation  and  an  equation  of  electron  current 
continuity,  each  with  specified  boundary  conditions. 

In  section  4.1,  the  statement  of  tnis  boundary  value  problem  is 
made  specific,  and  notation  is  presented.  The  employed  computation 
sch  me  is  describee  in  section  4.2,  and  a  semi -empi ri cal  algorithm  for 
automatically  calculating  a  suitable  sequence  of  relaxation  factors  is 
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given  In  section'  4.3.  The  difference  approximations  as  employed  are 
described  In  sections  4.4,  4.5;  some  additional  details  of  the  com¬ 
putations,  and  a  summary  of  numerical  experiments  designed  to  evaluate 
the  accuracy  of  the  computations,  are  given  In  section  4.6.  1 

Figure 63 Illustrates  the  two-dimensional  analytical  model  used 

i  i 

■  throughout  our  Investigation  of  bipolar  transistor  operation.  The 
boundary  conditions,  etc.,  used  to  approximate  this  structure  are 
discussed  in  section  4.7  and  the  associated  computational  method  is 
described  in  section  4.8. 

I 

i 

4.1  Definitions  and  Notation 

t  i 

We  employ  a  system  of  units  in  which  the  intrinsic  carrier  density, 
the  Boltzmann  voltage,  and  the  electronic  charge  have  magnitude  unity; 
lettingi  i[i(x,y) ,  n(x,y),  p(x,y)  denote,  respectively,  the  electrostatic 
potential,  electron  concentration,  and  hole  concentration  functions,  the 
Poisson  equation  becomes  i  .  i 

a)  esa^'-  Na  +  p  —  n  =  0,  (x,y)eR1 ;  ^ 

b)  Ai|i  3  0 ,  (x,y)eR2  ; 

where  es,  are  the  substrate  dielectric  constant  and  doping,1  respec¬ 
tively.  Across  the  substrate-oxide  interface,  ws  require  v  and  eaj'/By 

•  ,  * 

continuous;  the  other  boundary  conditions  are  shown  in  Fig.od.  Similarly, 

i  1 

the  equation  of 'electron  current  continuity  may  be  written 

vfpeMeA))3  0,  (x,y)eR1  (2) 


where  m  is  the  electron  mobility,  which  may  be  position  dependent  and 
anisotropic  due  to  dependence  on  the  local  electric  field.  On  the  sub¬ 
strate  boundary,  either  n  is  specified  (on  the  lines  CD  and  EF  in  Fig. 60) 

i 

or  the  normal  component  of  the  electron  current  density  vanishes.  In 
the  neglect  of  hole  current,  the  hole  concentration  can  be  written 


p  =  e 
where  Vgg 


is  the  source-substrate  or 


"back-gate"  bias  voltage. 


(3) 


In  the  following,  we  use  Vr,  Vg,  Vp,  respectively,,  fgr  the  built-in 
source  voltage,  the  applied  gate  voltage,  and  the  applied  drain  voltage. 
The  electron  quasi-Fermi  potential  is  set  to  zero  at  the  source. 
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We  Introduce  the  uniform  and  Dlrlchlet  norms  of  functions  defined 


In  R, 

a)  |  juj  j  *  Sup  |u(x,y)|  ; 
(x,y)sR 


4.2  Computation  Scheme 


In  this  section  we  describe  an  algorithm  for  computing  an  approxim¬ 
ation  to  a  simultaneous  solution  of  eqs.  1-3.  The  algorithm  Is  based 
on  an  Iteration  scheme  described  previously and  shown  to  converge 
globally  In  certain  cases.  Although  the  hypotheses  for  global  convergence 
are  not  satisfied  In  this  application,  the  Iteration  scheme  has  been 
found  entirely  satisfactory. 

We  first  compute  the  electrostatic  potential  on  the  line  ED  of 
Fig.  go  from  one-dimensional  theory 


S  3?^ 


na  + 


-*(x,K)-V, 


BG 


n(x,K) 


0,  xe(0,L) ; 


r  e^x,K),  x  <  L/2  , 

n(x,K)  =  \  *(x,K)-VD  *(5 

l  e  x  >  L/2  ; 

4.(0, K)  =  Vs,  *(L,K)  ■  Vs  +  VD  ; 

these  values  are  used  as  boundary  conditions  In  subsequent  calculations. 


An  Initial  approximation  to  the  electrostatic  potential  distri¬ 
bution  4.( 


Hx,y)  in  R  is  next  computed  from  the  Poisson  equation, 

,(0) 


-r 


=  0,  (x.yJeR- 


'BG  (0} 


n'  '(x,y)  =  0,  (>  ,y)eR  , 


(6) 


(o) 


together  with  the  appropriate  boundary  conditions.  In  (6),  n(x,y)  is 
an  electron  distribution  obtained  from  a  previous  computation  with  VQ 
and  Vgg  set  at  the  desired  values  and  Vp  at  a  lower  value  than  presently 
desired.  (If  necessary,  a  separate  computation  with  Vp  set  to  zero  is  made.) 
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The  following  Iteration  scheme  is  then  used  to  obtain  from  the 

(ml  Mj 

m-th  iterate,  <f>v  ;(x,y),  an  improved  approximation,  denoted  by  ip(x,y) 

and  the  corresponding  current  and  electron  distribution  n^(x,y). 

The  field  dependent  electron  mobilities  in  the  x-  and  y-directions, 

denoted  by  u  and  u  ,  are  calculated  from  the  equation  of 

A  Jr  t 

electron  current  continuity,  eq.  2,  Is  then  solved  for  n  '  by  des¬ 
cribing  the  current  density  components  in  terms  of  a  "stream  function" 
e(x,y)  defined  for  (x,y)eRj 


a) 


„  e**"0  i.  (e"*'(l">n(nl))  =  J(m)  ^ 
VXT  3X  n  '  u  3y 


(m) 


(7) 


b)  ,/(B)  <.-*W.M>  -  >>  ^ 

,(m) 


where  J'",/  is  the  m-th  approximation  to  the  total  device  current. 

Equations  7  may  be  readily  seen  to  imply  the  following  equation 
for  e( 


UUUU  I#  I 
(■) 

>(x,y) 


The  following  boundary  conditions  for  e^ni'  equivalent  to  the 
prescribed  boundary  conditions  for  n^,  in  view  of  the  relations  (7a,  7b) 

(m)  (m)  (m)  (m) 

9(x,0)  =  1;  e(x,K)  =  0;  f| - (y,0)  =  - (y,L)  -  0  .  (9) 

Equations  8,  9  are  solved  for  e|x,y);  the  value  of  is  then 
determined  by  integrating  (7a)  across  the  line  FC(y*0)  In  Fig. 60,' using 
the  source  and  drain  boundary  conditions  to  obtain 


e 


(m) 

-i|»(x,0)  de_ 

ay 


(m) 


<x,0) 


Then  the  electron  distribution  n 


(m) 


is  determined  from  (7a), 


(10) 


nlxiy)* 


(m) 

,-nz> 


%  *  ' 

b)  .e*^>  |«'V  /  «-*^»>  f-(z.y)  &]  . 

We  remark  that  this  method  of  solving  the  equation  of  current  con¬ 
tinuity  is  simply  a  generalization  of  the  one-dimensional  method  of 
Gummel;6  In  this  application,  however.  It  has  several  advantages:  the 
equation  for  the  current  distribution,  eq.  8,  can  be  solved  In  practice 
very  efficiently,  by  a  line  Iteration  described  below;  at  each  step, 
approximations  are  obtained  to  the  total  device  current  and  the  relative 
current  distribution;  the  use  of  appropriate  difference  approximations 
to  (11)  avoids  numerical  errors  due  to  roundoff  or  machine  overflow. 

At  each  step,  a  function  v^  defined  In  R  is  computed,  satisfying 
the  boundary  conditions  for  <^(x,y)  and  the  relations 


v e 


BS-  nW,  (x,y)cR,  ; 


-  0,  (xj')eR2  ;  ( 

and  the  Improved  approximation  to  f(x,y)  then  obtained  from 

(irH-1 )  (m)  (m)  (m)  (m)  (m-1) 

</'(x,y)=v(x,y)+am(v{x,y)-^(x,y))+3m(i/<(x,y)-i(;(x,y)),  (x,y)eR  ( 

where  a,  are  constants  evaluated  below, 
m  m 

The  Iteration  is  terminated  when  the  valec  of  |  |v^-<j/rn^  1 1  falls 
below  a  specified  value. 

4.3  Relaxation  Factors 

Suitable  expressions  for  the  constants  am,  3m  appearing  in  (13) 
are  obtained  by  making  the  same  linearization  assumption  originally 
described  by  Gummel,6  and  using  the  Chebyshev  sequence  of  relaxation 

7 

factors . 
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Let  u(m>  =  4)  -  be  the  error  in  the  m-th  iterate,  and  let 
S  denote  the  set  of  functions  which  are  twice  continuously  differen¬ 
tiable  in  Rp  harmonic  in  Rg,  vanish  on  the  boundary  of  R  (the  lines 
ABDEA  of  Fig.  1),  and  satisfy  the  interface  relations  prescribed  for 
<p(x,y)>  i.e.,  u(x,0+)su(x,0-) ,  es~(x,G+)=eox|y<x,0-),  for  a11  UtS- 

Clearly,  u^eS.  We  now  assume  u^  is  small,  and  that  the 
dominant  mechanism  affecting  the  propagation  of  the  error  u^  is 

fi  R 

the  exponential  dependence  of  the  carrier  densities  on  the  local * 
electrostatic  potential,  as  shown  in  eqs.  3,  11.  Then  from  the  mean 
value  theorem,  we  obtain  the  approximation 

n(m)  _  p(m)  -  n  _  p  .  (n+p)  u(m)f  (14) 

where  n,p  denote  the  electron  and  hole  concentrations  in  the  desired 
solution. 

It  is  convenient  to  expand  u^  in  the  set  of  eigenfunctions  Zj, 
corresponding  to  eigenvalues  A.  satisfying 

J 

xjesAzj  +  ^n+p^  2W  =  °*  (x»y)eRi  »  zjeS  •  (is) 

The  z,  form  a  complete  set  in  S  and  satisfy  the  following 
J  q 

orthogonality  relation 


/ (n+p)  Zj(x.y)  zk(x,y)  dx  dy  ■  «.  k. 


(16) 


Subtracting  eq.  la  from  eq.  12  gives 
e  a(v^-i|>)  =  n^-  p^-  n  +  p  =  -(n+p)u^  , 


07) 


using  the  approximation  14.  It  follows  from  15,  16,  17  that  the  'z^ 
are  also  eigenfunctions  of  the  mapping  u^+v^-/11^,  with  corresponding 
eigenvalues  -(1+Aj).  Denoting  by  a  a  bound  on  the  Aj ,  which  are 
positive  and  accumulate  only  at  zero,  it  follows  (see  ref.  7  for  details) 
that  appropriate  values  of  am  and  Bm  are  given  by 


a) 


_  4  cosh  kw 
“m"  A  cosh  (k+l )u>  *  e 


cosh  (k-1 
m  cosh 


Ms ;  c°sh 


=  1  +  2/A  ,  k  f  0; 


08) 


b)  =  2/a,  =  0,  k  =  0. 

In  18  a,  k  is  a  nonzero  integer  depending  on  the  iteration  number  m. 
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In  each  computation,  k  Is  set  to  zero  Initially  and  Increased  by 
one  at  each  Iteration.  However,  at  any  Iteration  for  which 

>  1  j  { ,  the  value  of  k  is  reset  to  zero, 

such  a  condition  being  considered  an  Indication  that  nonlinear  effects 
have  become  important. 

It  remains  to  estimate  the  constant  X;  It  follows  from  eq.  15 
that  x  Is  given  by 


J"  (n+p)z^dx  dy 


x  =  Sup 
ZeS 


The  right  side  of  (19)  may  be  estimated  by  several  methods.  For 
example.  It  follows  from  Holder's  Inequality  and  the  Sobolev  lemma 
that  for  any  q  >  1 ,  there  exists  a  constant  C(q)  such  that 


a  i  c(d)  f  (n+p)q  dx  dy 


However,  the  sharply  peaked  form  of  the  electron  distribution 
function  at  the  oxide  interface  makes  estimates  of  this  form  too  con¬ 
servative  to  be  of  much  practical  value.  Although  eq.  20  is  not 
strictly  valid  for  q  =  1,  we  have  found  empirically  that  a  useful 
value  of  A  is  given  by 


*  ~  4e  L  f  (n+p)  dxdy  . 


In  evaluating  eq.  21,  it  is  also  necessary  to  replace  the  unknown 
functions  n,p  with  the  initial  approximations  n^,  p^°\ 

4.4  Difference  Approximations 

The  computations  described  above  are  readily  performed  on  a  computer 
by  the  use  of  finite  difference  approximations. 
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In  view  of  the  discontinuity  of  the  dielectric  constant  at  the 
oxide  interface  *  it  is  convenient  to  discretize  the  Tolsson  equation 
in  terms  of  the  variables  x,z,  where 


z=z(y)= 


y/'ox*  y  i  °* 
y/es .  y  1  o  ; 


(22) 


the  Poisson  equation  then  becomes 


^f  +  (NA+  P'n>’  2  ’  °  ! 


(23) 


since  and  ||-  are  continuous  across  the  oxide  interface  (z=0),  no 
additional  relations  there  are  necessary. 


We  superimpose  a  rectangular  grid  on  the  geometry  of  Fig. 60;  the 
locations  of  the  grid  lines  x^,  yj  are  controlled  by  the  following 
relations: 


a) 

b) 

c) 


xi+l  "  X1  5  hx(1+1)  =  X)hx(i) 


yj+l  ’  yj  H  hy^  =  e  hz(j)  = 


(l'6y)hy(j-l) 

(HfiyJhyU-!) 


hz(m)  «  (l-6y)nz(m-l)  =  hz(m+l)/(l+6y) 


j  <  m 
j  >  nH-1 


(24) 


d)  hx(i)  =  l/2(hx(1)+hx(1+l),  hy(j)  =  l/2(hy ( j-1 )+hy(j ) ) 

where  m  is  the  y-index  corresponding  to  the  oxide  interface  and  the 

constants  5  ,  6  are  defined  below, 
x  y 

The  orientation  of  the  discrete  variables,  relative  to  the  grid 
points,  is  shown  in  Fig.  61,  and  the  treatment  of  the  boundary  conditions 
in  Figs. 62a  and 62).  The  orientation  of  the  discrete  variables  near  the 
oxide  interface  is  shown  in  Fig. 62b.  drawn  approximately  to  scale  for 
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1 

0  (l,j) 

1 
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1 

f 
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l/2[hy(H)-<hy<j)]  j  l  | 

^  My(i-ij-i)  e a, j-ii)  My(ij-i)  ,  ea^j-D  My<i+ij-i) 
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n(i-l,j-l)  1  n(i,  j-1)  n(i+l,j-l) 
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•*hx(i!=  2  Ihx(1)+hx 11+1)1  -H 


Fig.  61.  Orientation  of  discrete  variables. 
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Fi«g.  62a.  Boundary  conditions  at  source  contact. 
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Fig.  62b.  Boundary  conditions  at  oxide  interface. 


By  the  standard  "Une  Integral"  method,  a  difference  approximation 

to  eq.  23  Is  obtained  whose  local  truncation  ^ror  Is  of  the  order 
2  2 

h„  +  6vhv  +  h,  +  he.  By  the  same  method,  we  obtain  difference 
a  a  a  z  z  y 

approximations  to  the  other  required  equations  (8,  10,  11}.  Using 
a(i+l/2,j),  a(i,j+l/2),  respectively,  as  values  of  between  the 
pairs  of  points  (i,jj  1+1  ,j)  and  (i,j;  i,j+l),  we  obtain  the  folding; 


where  the  iteration  denoting  superscripts  and  the  difference  equation 
approximating  to  the  Poisson  equation  have  been  omitted  for  simplicity. 


The  systems  of  linear  equations  generated  by  these  difference  ap¬ 
proximations  can  all  be  efficiently  solved.  The  coefficient  matrix 
corresponding  to  the  Poisson  difference  equations  is  symmetric  and  ; 

positive  definite,  and  can  readily  be  solved  by  the  method  of  Dupont,  t 

10  1 
Kendall,  and  Rachford.  To  solve  the  system  generated  by  eq.  25,  we  ; 

exploit  the  fact  that  the  current  flow  in  the  I6FET  device  is  essentially  > 

one-dimensional,  i.e.  the  horizontal  gradient  terms  in  eq.  25  are  I 

11  j 

relatively  unimportant.  We  thus  employ  a  line  iteration,  solving  ! 

successive  tri diagonal  systems  of  m  equations  generated  by  eq.  25  for  j 

a  fixed  i,  and  j  varying.  For  this  particular  model,  this  iteration  j 

converges  very  rapidly.  Equation  26  is  explicit  for  J,  and  27  for  n;  j 

in  practice,  we  solve  eq.  27  for  n(i+l,j)  using  the  value  of  n(i,j)  if  j 

^(i+l,j)  <  ^(i,j),  and  vice  versa.  This  procedure  suppresses  the  buildup 
of  roundoff  error;  the  consistency  of  the  difference  equations  (25,26,27)  * 

insures  the  continuity  of  t*»e  obtained  values  of  n(i,j). 


Ill 


4.5  A  Critical  Factor 


It  remains  to  define  the  quantities  a(1+l/2,j),  a(1,j+l/2)  in 
terms  of  local  values  of  the  electrostatic  potential;  we  employ  the 
following  approximation; 


e-*(1+U)  _e-*(1J) 


a(i+l/2,j)~ 


,j)  v*(v+r,'ji -  * 


(20) 


-  *0+1 J). 


with  an  analogous  expressing  for  a( 1 , j+1/2) .  This  expression  m«y  be 

12 

derived,  as  was  done  by  Scharfetter  and  Gunnel ,  by  noting  that 
a(i+l,j)  is  identified  with  the  average  value  of  the  function  exp(-i/0 
in  the  interval  xe(x.j,  x^),  y  *  y^,  and  assuming  a  linear  dependence 
of  t|»  on  x  in  this  interval.  This  association  becomes  apparent  in  eq.  27, 
if  one  anticipates  variations  of  several  units  of  kT/q  between  i|»(1,j) 
and  *(i+1  J)»  the  left  side  of  eq.  27  then  cannot  be  regarded  as  an 
accurate  difference  approximation  to  hx  —  (ne~^)  in  this  Interval. 

One  is  motivated  to  Interpret  eq.  27  as  an  Integrated  form  of  7a  bet¬ 
ween  the  points  (x^,  y^)  and(x^+-j,  yj). 

In  one -dimensional  models,  it  can  also  be  shown  by  purely  analytic 
13 

arguments  that  the  expression  (23)  leads  to  a  discretization  error 
of  0(h)  in  the  computed  device  current,  without  assuming  small  variations 
in  ip  between  successive  grid  points. 


Finally,  we  show  that  an  obvious  alternative  to  eq.  28,  which  has 

2  3 

been  used  in  previous  calculations  on  IGFET  models, 

a*(i+l/2,j)  =  exp  j  (i|»(i,j)  +  *(  i  «•! ,  j )  )j  ,  (29) 

can  lead  to  serious  error  in  the  computed  electron  distribution. 

Suppose  (x.j,  yj)  is  u  point  near  the  drain  junction  at  which  a 
significant  current  density  exists,  in  a  device  operating  in  its 
saturation  region.  We  assume  that 

*0+1  J)  -  *OJ)  »  1  (30) 


and  consider  the  solution  of  eq.  27  for  n(i,j).  Setting 


J(i+l/2,j)=J  iJ."1)  ,  E(i+l/2,j)  = 

y  ^ 

the  local  current  density  and  electric  field  In  the  x-direction,  res¬ 
pectively,  we  obtain  from  eqs.  27,  28,  30, 

n(1,j)  *  y  (l+(  J)ffi+jy2Vj)  *  0texP^^i »J)-*C1+1  »3)))i  (31) 

the  '•cmputed  electron  density  Is  sufficient  to  support  the  local  current 
density. 

We  can  estimate  the  direction  of  the  error  in  eq.  31.  From  the 

2  2 

Poisson  equation,  it  is  clear  that  3  <|>/ sx  is  positive  in  such  an 
interval  x^  <  x  <  x^+, ;  therefore  the  function  e"^  is  larger  in  this 
interval  than  the  linear  approximation  gives,  and  our  expression  for 
a(i+l/2,j)  is  an  underestimate.  It  follows  that  eq.  31  is  an  under¬ 
estimate  of  the  electron  concentration. 

If,  however,  eq.  29  is  used  instead  of  eq.  28,  we  obtain  an 
electron  concentration  n*(i,j)  given  by 


n*(1J)  «  n(i,j) 


feu 

l+Lll 

l-<|/(i,j) 

) 

l  i 

sinh 

»(1+1,J^-»(i 

31) 

(32) 


where  n(i,j)  is  the  electron  concentration  obtained  in  eq.  31.  The 

value  of  n*(i,j)  can  thus  be  unrealistically  small;  for 

</'(i+1  .j)-ij>(1  »j)  -  20  (0.5  volt),  we  obtain  n(i ,j)/n*(i ,j)  =  103. 


4.6  Summary  -  Numerical  Experiments 

A  unified  presentation  of  the  IGFET  device  characteristics  as 
computed  by  the  method  described  above,  including  a  detailed  comparison 
with  existing  theory  appears  elsewhere  in  this  report.  The  empirical 
data  given  here  is  intended  to  support  the  numerical  accuracy  of  the 
results  presented. 


Most  of  our  computations  have  been  run  with  20  horizontal  and 

30  vertical  grid  points,  with  6  =  0.05,  6  -  0.1  in  eq,  24,  Empirical 

x  y 

accuracy  estimates  have  been  obtained  by  comparison  of  the  results 
with  results  obtained  using  higher  point  densities  and  different  values 
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of  6X,  6^;  on  this  basis,  the  numerical  (discretization)  errors  in 
the  presented  results  are  estimated  at  5*  of  the  quoted  values  of 
device  currents,  and  1$  +  20  mV.  of  the  quoted  values  of  applied  , 
voltage  along  the  substrate>ox1de  Interface.  , 

Numerical  experiments  have  also  been  conducted  to  assess  the  , 
Importance  of  the  difference  approximation  (2(5),  the  reliability  of 
the  formula  (21),  and  the  sensitivity  of  the> computed  results  to  the 
precise  condition  adopted  for  terminating  a  computation.  The  predlc- 
tlons  of  section  4.5  regarding  the  computed  electron  distribution  have 
been  observed  in  practice. 

I 

The  use  of  a  Chebyshev  sequence  cf  relaxation  factors,  as  given 
by  eq.  18,  reduces  the  sensitivity  of  the  number  of  Iterations  required 
in  a  computation  to  the  sharpness  of  the  adopted  value  of  the  constant 
A  and  the  stopping  condition.  The  relation  (21)  has  been  found  suf¬ 
ficiently  conservative  to  insure  convergence  In  all  cases  tested;  in 
a  few  test  cases.  It  was  found  to  give  a  value  of  a  wlthlh  a  factor  of  1 
two  of  the  optimum.  > 

The  relaxed  Picard  Iteration  scheme  described  above  converges 
essentially  in  the  norm  1 1 1 1 1  ;5  previous  analysis^3  suggests 
that  in  this  case,  the  computed  value  of  the  device  current  will  not 
depend  strongly  on  the  adopted  terminating  condition.  Our  experiments 
confirmed  this,  and  the  results  given  in  this  report  were  run  using 

||v(m)  -  <  10  mV.  11  •  (33) 

as  the  stopping  criterion.  We  remark  that  this  property  is  particular 
to  the  Picard  iteration  scheme.  i 

4.7  Bipolar  Transistor  Model  1 

,  i 

A  numerical  analysis  of  the  bipolar  transistor  has  also  been  ; 
developed  and  incorporated  into  a  computer  program.  In  this  section 
we  present  the  associated  equations  and  boundary  conditions;  the 
computation  scheme  is  described  in  section  4.8. 

The  geometry  associated  with  our  approximation  for  the  bipolar 
transistor  is  shown  in  Fig. 63.  Letting  n  denote  the  rectangle 
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Emitter  Contact 


K' 


Fig.  63.  Analytical  model  for  a  bipolar  transistor 


I  I 

o<x<K,o<y<L,  we  consider  the  electrostatic  potential  ^  and 
the  local  carrier  densities  n,  p,  to  be  described  by  the  Poisson  equation 
and  continuity  equations, 

eAt|>  +  N(x,y)  +  p  -  n  *  0,  (34) 

v  =  Rn(n,p)»  (35) 

v  (jipe“*7(fc*p))  =  Rp(n.p);  (x,y)  in  q  .  (36) 


In  eq.  34,  e  Is  the  dielectric  constant  and  N(x,y)  the  specified 
doping  concentration,  which  determines  the  junction  locations.  In 
these  calculations,  we  have  used  profiles  of  the  form 

-kv2 

N(x,y)  -  C£f(x,y)  -  Cbe  y  +  Cc  ,  (37) 

where  Ce,  Cb,  Cc,  k  are  constants  and  f(x,y)  Is  the  solution  of  a 
two-dimensional  emitter  diffusion  continuity  equation.  In  the  case 
of  a  constant  Impurity  diffusion  coefficient,  a  series  solution  to 
this  problem  Is  known;14  alternatively,  the  numerical  integration 
of  this  equation  Is  not  difficult.  For  simplicity,  we  have  assumed 
In  these  calculations  that  the  emitter  contact  and  emitter  diffusion 
mask  opening  coincide. 


In  eqs.  35,  36,  we  have  assumed  nondegenerate  statistics;  we 
also  have  used  the  standard  Hall -Shockley-Read  recombination  term, 

-  1 


Rn(n’P}  “  Rp(n,p)  =  Tn(n+??+%pIp^) 


(38) 


The  electron  and  hole  mobilities  in  eqs.  35,  36  are  dependent 
on  position,  local  electric  field,  and  direction.  Position  dependency 
is  Included  by  the  method  of  Caughey  and  Thomas.15’  In  the  collector 
junction  space-charge  region,  where  current  flow  is  assisted  by  the 
local  electric  field,  a  correction  to  the  carrier  mobilities  is  applied 
by  relations  of  the  form 


PnjX(x,y)  =  minimum 


Vo** •*>. 


where  un  (x,y)  is  the  x-directiori  electron  mobility  at  the  point 

li  jA 

(x,y),  un  0(x,y)  is  the  position-dependent  electron  mobility  at  (x,y). 
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E  Cx *y)  the  local  x-dtrectton  component  of  the  electric  field,  and  Vr  _, 
Vj  n  are  the  "critical"  and  "limiting"  electron  velocities  respectively. 
Similar  relations  hold  for  un^y(x,y),  up>x(x,y),  pp^y(x,y),  with  the 
obvious  changes  in  subscripts. 

In  the  emitter  junction  space-charge  region,  current  flow  is 
opposed  to  the  local  electric  field,  and  the  mechanisms  of  carrier 
transport  are  not  well  understood.  In  this  region,  we  have  used 
relations  of  the  form  (eq.  39)  with  the  "critical  velocity"  term 
omitted. 

Boundary  conditions  for  the  system  (eqs.  34-36)  are  formulated  as 
follows:  the  contact  boundary  segments,  AF,  BC,  OE  in  Fig.  63,  are 
assumed  to  be  infinite  surface  recombination  velocity  ohmic  contacts, 
at  which  thermal  equilibrium  and  charge  neutrality  apply;  thus  n,  p, 
and  t|>  can  all  be  specified.  The  applied  bias  voltages  ere  defined  in 
the  usual  way,  in  terms  of  the  carrier  quasi -Fermi  potentials.  Along 
the  boundary  segments  AB,  CD,  EF,  the  normal  component  of  the  electron 
and  hole  current  density  is  required  to  vanish.  On  the  boundary  seg¬ 
ments  AB  and  EF,  the  electrostatic  potential  distribution  is  computed 
from  one-dimensional  theory  and  specified,  analogously  to  the  I6FET 
model,  eq.  5;  on  the  segment  CD,  the  normal  component  of  the  electric 
field  is  also  required  to  vanish. 

4.8  Computation  Scheme  for  Bipolar  Transistor  Model 

The  computation  scheme  described  in  section  4.2  is  inappropriate 
for  the  bipolar  transistor  model  for  two  reasons:  the  finite  recom¬ 
bination  term  precludes  the  stream  function  treatment  of  the  carrier 
continuity  equations,  and  the  large  carrier  densities  would  reduce 
the  admissible  relaxation  factors  to  very  small  values,  as  shown  in 
eq.  21.  We  have  adopted  the  following  variation  of  Gunnel's 
algorithm.6 

At  each  point  (x,y)  in  n,  we  assume  that  approximations  to  the 
electrostatic  potential  and  the  recombination  term  are  known,  denoted 
by  </m^(x,y),  R^(x,y)  respectively.  We  compute  corresponding 
approximations  to  the  carrier  densities  and  "improved"  approximations 


to  the  potential  and  recombination  terms  by  successively  solving  the 
following  equations,  each  with  its  appropriate  boundary  conditions, 
as  described  in  section  4.7. 


*(m))  -  e(n(m)+  plm))U(inH)-  *(m)) 

+  a(£4*^+  N(x,y)  +  n*mh  *  0  B3) 

In  eq.  43,  a,  e  are  positive  constants  of  order  unity,  suitable 
values  of  which  have  been  determined  empirically  for  each  computation 
reported.® 

Difference  approximations  to  eqs.  40-43  are  readily  obtained, 
by  methods  similar  to  those  used  for  the  IGFET  analysis,  and  described 
in  section  4.4.  The  exponential  average  factors,  described  in  section 
4.5,  have  also  been  incorporated  in  the  bipolar  device  computations. 

The  computations  described  in  this  report  were  run  with  40 
horizontal  and  70  vertical  grid  points.  The  numerical  accuracy  of  the 
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scheme  was  appraised  by  rerunning  we  computation,  using  a  75  x  120 
point  grid.  The  two  values  of  collector  current  and  hpE  obtained  in 
this  manner  varied  by  2%  and  4%,  respectively.  Typical  variation 
In  the  local  carrier  densities  v/as  2X,'and  in  the  electrostatic 
potential  distribution  was  5  mV.  We  therefore  anticipate  numerical 
errors  of  this  magnitude  in  the  results  presented. 
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